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Pion, Kaon, Proton and Antiproton Spectra in d+Au and p+p Collisions 
at ^/sNN= 200GeV at the Relativistic Heavy Ion Collider 
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Jan. 2005 

Identified mid-rapidity particle spectra of vr^, K^, and p{p) from 200 GeV p+p and 
d+Au collisions are reported. The d+Au collisions were divided into 3 centralities. 
This data were taken in 2003 run from the Solenoidal Tracker at RHIC (STAR) 
experiment. A time-of-flight detector based on multi-gap resistive plate chamber 
(MRPC) technology is used for particle identification. This is the first time that 
MRPC detector was installed to take data as a time-of-flight detector in the collider 
experiment. 

The calibration method was set up in the STAR experiment for the first time, 
which will be described in the analysis chapter in detail and has been applied to the 
experimental data successfully. The intrinsic timing resolution 85 ps was achieved 
after the calibration. In 2003 run, the pion/kaon can be separated up to transverse 
momentum (px) 1.6 GeV/c while proton can be identified up to 3.0 GeV/c. Thus the 
identified particle spectra can be extended to intermediate pt in STAR. This proved 
that MRPC as a time-of-flight detector works in the heavy ion collider experiment. 

We observe that the spectra of vr^, K^, p and p are considerably harder in d+Au 
than those in p+p collisions. In y/s^ = 200 GeV d+Au collisions, the nuclear 
modification factor RdAu of protons rise faster than those of pions and kaons. The 
RdAu of proton is larger than 1 at intermediate pt while the proton production follows 
binary scaling at the same pt range in 200 GeV Au+Au collisions. These results 
further prove that the suppression observed in Au+Au collisions at intermediate and 
high Pt is due to final state interactions in a dense and dissipative medium produced 



during the collision and not due to the initial state wave function of the Au nucleus. 
Since the initial state in d+Au collisions is similar to that in Au+Au collisions, and, 
it's believed that the quark-gluon plasma doesn't exist in d+Au collisions, these 
results from d+Au collisions are very important for us to judge whether the quark- 
gluon plasma exists in Au+Au collisions or not and to understand the property of the 
dense matter created in Au+Au collisions. Besides, the particle-species dependence 
of the Cronin effect is observed to be significantly smaller than that at lower energies. 
The ratio of the nuclear modification factor (RdAu) between {p + p) and charged 
hadrons (h) in the transverse momentum range 1.2 < pt < 3.0 GeV/c is measured to 
be 1.19 ± 0.05(stat)±0.03(syst) in minimum-bias collisions and shows little centrality 
dependence. The yield ratio of {p + p)/h in minimum-bias d+Au collisions is found 
to be a factor of 2 lower than that in Au+Au collisions, indicating that the relative 
baryon enhancement observed in heavy ion collisions at RHIC is due to the final state 
effects in Au+Au collisions. 

The mechanism for Cronin effect is also discussed in this thesis by comparison 
with the recombination model and with the initial multiple parton scattering 
model IHO], which will be described in the discussion chapter in detail. Usually the 
Cronin effect has been explained to be the initial state effect only [30j. However, 
from the comparisons, we conclude that the Cronin effect in ^,/s^ = 200 GeV d+Au 
collisions is not initial state effect only, and that final state effect plays an important 
role. 

These physics results has been at e-Print Archives (nu-ex/0309012) and submitted 
for publication. The excellent particle identification from the prototype MRPC tray 
and the important physics from it have provided a solid basis for the STAR full-time- 
of-flight-system proposal. 
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Chapter 1 
Physics 

1.1 Deconfinement and phase diagram 

The theory which describes the interaction of the color charges of quarks and glu- 
ons is called Quantum Chromodynamics (QCD). In phenomenological quark models, 
mesons can be described as quark- ant iquark bound states, while baryons can be con- 
sidered as three quark bound states. Up to now, it's found that all the hadron states 
which can be observed in isolation is colorless singlet states. Experimentally, no sin- 
gle quark, which is described by a color-triplet state, has ever been isolated. The 
absence of the observation of a single quark in isolation suggests that the interac- 
tion between quarks and gluons must be strong on large distance scale. In the other 
extreme, much insight into the nature of the interaction between quarks and gluons 
on short distance scales was provides by deep inelastic scattering experiments. In 
these experiments, the incident electron interacts with a quark within a hadron and 
is accompanied by the momentum transfer from the electron to the quark. The mea- 
surement of the electron momentum before and after the interaction allows a probe of 
the momentum distribution of the parton inside the nucleon. It was found that with 
very large momentum transfer, the quarks inside the hadron behave as if they were 
almost free jT]. The strong coupling between quarks and gluons at large distances 
and asymptotic freedom are the two remarkable features of QCD. When the energy 
density is high enough either due to the high temperature or high baryon density. 



1 



2 



the quark or gluon may be deconfined from a hadron. The thermahzed quark gluon 
system is what we called quark-gluon plasma. Lattice QCD calculations, considering 
two light quark flavors, predict a phase transition from a confined phase, hadronic 
matter, to a deconfined phase, or quark-gluon plasma (QGP), at a temperature of 
approximately 160 MeV Figure 1.1 shows the phase diagram of the hadronic and 

I ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' 

early universe 
I 

LHC 

I quark-gluon 
RHIC plasma 




baryonic chemical potential [GeV] 



Figure 1.1: Phase diagram of hadronic and partonic matter. Figure is taken from 0. 

partonic matter. A phase transition from the confined hadronic matter to the decon- 
fined QGP matter is expected to happen at either high temperature or large baryon 
chemical potential fiB- Recent Lattice QCD calculations show that the QGP is far 
from ideal below 3 Tc. The nonideal nature of this strongly coupled QGP is also seen 
from the deviation of the pressure, -P(T), and energy density e(T) from the Stefan 
Boltzmann limit as shown in Figure from Experiments on relativistic heavy ion 
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1 1.5 2 2.5 3 1 1.5 2 2.5 3 

T/T, T/T, 

Figure 1.2: A recent Lattice QCD calculation of the pressure, P(T)/T^, and a 
measure of the deviation from the ideal Stefan-Boltzmann limit (e(T) — 3P{T))/T^. 

collisions are designed to search for and study the deconfined QGP matter. 

1.2 Relativist ic Heavy Ion Collisions 

The experimental programs in relativistic heavy ions started in 1986 using the Al- 
ternating Gradient Synchrotron (AGS) at Brookhaven National Lab (BNL) and the 
Super Proton Synchrotron (SPS) at European laboratory for particle physics (CERN). 
At BNL, ion beams of silicon and gold, accelerated to momenta of 14 and 11 GeV/c 
per nucleon, respectively, have been utilized in 10 fixed-target experiments. There 
have been 15 heavy ion experiments at CERN utilizing beams of oxygen at 60 and 
200 GeV/c per nucleon, sulphur at 200 GeV/c per nucleon and Pb at 160 GeV/c per 
nucleon j2]. 

The Relativistic Heavy Ion Collider (RHIC) at BNL is designed for head-on Au+Au 
collisions at ^/'sNN = 200 GeV. The first RHIC run was performed in 2000 with 
Au+Au collisions at y/sj^ = 130 GeV/c in four experiments, STAR, PHENIX, PHO- 
BOS and BRAHMS. The second RHIC run was in 2001 and 2002 with Au+Au and 
p+p collisions at ^/sNN = 200 GeV. The third RHIC run was in 2002 and 2003 with 



4 



d+Au and p+p collisions at ^/snn — 200 GeV. 

The above mentioned relativistic heavy ion collision experiments are designed for the 
search and study of the possible deconfined high energy density matter, quark-gluon 
plasma. In head-on relativistic heavy ion collisions, two nuclei can be represented as 
two thin disks approaching each other at high speed because of the Lorentz contrac- 
tion effect in the moving direction. During the initial stage of the collisions, the energy 
density is higher than the critical energy density from the Lattice QCD calculation, 
so the quarks and gluons will be de-confined from nucleons and form the quarks and 
gluons system. The large cross section of interaction may lead to the thermalization 
of the quarks and gluons system. That's what we called the formation of quark-gluon 
plasma (QGP). In this stage, the high transverse momentum (p^) jets and cc pair 
will be produced due to the large momentum transfer. After that, the QGP will ex- 
pand and cool down and enter into the mixed-phase expansion. The chemical freeze 
out point will be formed after the inelastic interactions stop. That means that the 
particle yields and ratios will not change. After the chemical freeze out, the elastic 
interactions between hadrons will change the pt distribution of particles. The parti- 
cles will freeze out finally from the system after the clastic interactions stop. That's 
what we called the kinetic freeze out point. In the following the important results 
from RHIC will be addressed. 



1.3 The experimental results at RHIC 
1.3.1 Flow 

In non-central Au+Au collisions, the spatial space asymmetry will be transferred into 
the momentum space asymmetry by the azimuthal asymmetry of pressure gradients. 
The azimuthal particle distributions in momentum space can be expanded in a form 
of Fourier series 

n=l 
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Transverse Momentum pj (GeV/c) 

Figure 1.3: The minimum-bias (0-80% of the coUision cross section) V2{pt) for Kg, 
A + Aand/t±. The error bars shown include statistical and point-to-point systematic 
uncertainties from the background. The additional non-flow systematic uncertainties 
are approximately -20%. Hydrodynamical calculations of f 2 for pions, kaons, protons 
and lambdas are also plotted Figure is taken from [H]. 

where denotes the reaction plane angle. The Fourier expansion coefficient Vn 
stands for the nth harmonic of the event azimuthal anisotropies. vi is so called direct 
flow and f 2 is the elliptic flow. The elliptic flow is generated mainly during the highest 
density phase of the evolution before the initial geometry asymmetry of the plasma 
disappears. Hydrodynamical calculations [IHI show most of V2 is produced before 3 
fm/c at RHIC. 

Figure IT7S1 shows that The t>2 of Kg, A + A and charged hadrons (h^) as a function of 
Pt for 0-80% of the collision cross section ,6J. Also shown are the V2 of pions, kaons, 
protons and lambdas from hydrodynamical model The V2 from hydrodynamical 
model shows strong mass dependence, which fits the Kg V2 up to ~ 1 GeV/c and 
fits the A + A f 2 up to ~ 2.5 GeV/c. Even though the V2 from hydrodynamical 
model shows consistency with data at low Pt, however, the V2 from experimental 
results show saturation at intermediate px while hydrodynamical predictions show 
rising trend at the same pt range. 



6 



1.3.2 High pt suppression and di-hadron azimuthal correla- 
tion 

The V2 from hydrodynamical models show consistency with data at lower px and 
fail to reproduce data at higher pt- At high pt, the suppression for charged hadron 
production was observed in Au+Au collisions at RHIC energy. The comparison of the 
spectra in Au+Au collisions through those in p+p collisions, scaled by the number 
of binary nucleon nucleon collisions is the nuclear modification factor Raa- 

where Taa = (Abin)/o'hiii' accounts for the collision geometry, averaged over the event 
centrality class. ( A^bin) , the equivalent number of binary A^A^ collisions, is calculated 
using a Glauber model. The Raa is an experimental variable. The high pj- hadron 
suppression in central Au+Au collisions can also be investigated by comparing the 
hadron spectra in central and peripheral Au+Au collisions. That's what we called 
Rcp- Rcp is defined as 

{N^in^'^^)d'^NP''''P^''''^ydpTdr]' ^ ' ' 



Figure ^3 shows Raa{pt) of inclusive charged hadron for various centrality bins 
in Au+Au collisions at ^JsMN=2QQ GeV. Raa{pt) increases monotonically for px < 
2 GeV/c at all centralities and saturates near unity for p^ > 2 GeV/c in the most 
peripheral bins. In contrast, Raa{pt) for the central bins reaches a maximum and 
then decreases strongly above px = 2 GeV/c, showing the suppression of the charged 
hadron yield relative the A^A^ reference jH]. 

Suppression of high px hadron production in central Au+Au collisions relative to 
p+p collisions [HI E] has been interpreted as energy loss of the energetic partons 
traversing the produced hot and dense medium 0, that's so called jet quenching. If a 
dense partonic matter is formed during the initial stage of a heavy-ion collision with 
a large volume and a long life time (relative to the confinement scale I/Aqcd), the 
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Figure 1.4: Raa{pt) of inclusive charged hadron for various centrality bins. Figure is 
taken from 0. 



produced large Et parton will interact with this dense medium and will lose its energy 
via induced radiation. The energy loss depends on the parton density of the medium. 
Therefore, the study of parton energy loss can shed light on the properties of the 
dense matter in the early stage of heavy- ion collisions 9J. At sufficiently high beam 
energy, gluon saturation is also expected to result in a relative suppression of hadron 
yield at high px in A+A collisions [33 . Also shown in the Figure 11.41 are the results 
from perturbative QCD (pQCD) calculations. The Full-pQCD calculations include 
the partonic energy loss, the Cronin enhancement (due to initial multiple scattering) 
and the nuclear shadowing effect. The suppression is predicted to be pt independent 
when pt is larger than 6 GeV/c, which is consistent with our data. However, the 
discrepancy at 2-6 GeV/c was observed between the prediction and the experimental 
data. This discrepancy may be due to different mechanism for particle production at 
intermediate pt- The particle production at intermediate pt will be discussed later 
in this chapter. 
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Figure 1.5: (a) Azimuthal distribution of particles with respect to a trigger particle for 
p+p collisions (solid line), and mid-central Au+Au collisions within the reaction plane 
(squares) and out-of-plane (circles) at 200 GeV ^Uj. (b) Mean transverse momentum 
for particles around the away-side region as a function of number of charged parti- 
cles iHjj. The solid line shows the mean transverse momentum of inclusive hadrons. 



A more differential probe of parton energy loss is the measurement of high px di- 
hadron azimuthal correlation relative to the reaction plane orientation. The trigger 
hadron is in the range 4 < pt < Q GeV/c and the associated particle is at 2 < < 4 
GeV/c. Figure IT7S1 (left) shows the high pt di-hadron correlation when the trigger 
particle is selected in the azimuthal quadrants centered either in the reaction plane 
(in plane) or orthogonal to it (out of plane). The near side di-hadron azimuthal cor- 
relations in both cases were observed to be the same as that in p+p collisions, while 
the suppression of back to back correlation shows strong dependence on the relative 
angle between the triggered high px hadron and the reaction plane. This systematic 
dependence is consistent with the picture of parton energy loss: the path length for 
a dijet oriented out of plane is longer than that for a dijet oriented in plane, leading 
to a stronger suppression of parton energy loss in the out of plane. The dependence 
of parton energy loss on the path length is predicted to be substantially larger than 
linear jH]. 

The energy lost by away side partons traversing the collision matter must in the form 
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of the excess of softer emerging particles due to the transverse momentum conserva- 
tion. An analysis of azimuthal correlations between soft and hard particles has been 
performed for both 200 GeV p+p and Au+Au collisions [TJ at STAR as a first of 
attempt to trace the degree of the degradation on the away side. With triggered 
hadron still in the range 4 < p^p^ < 6 GeV/c, but the associated hadrons now sought 
over 0.15 < p-r < 4 GeV/c, combinatorial coincidences dominate this correlation and 
they must be subtracted carefully by mixed-event technique and also the elliptic flow 
effect was also subtracted by hand The results demonstrate that, in comparison 
with the p+p and peripheral Au+Au collisions, the momentum-balancing hadrons 
opposite to the high pt triggered particle in central Au+Au are greater in number, 
much more widely dispersed in azimuthal angle, and significantly softer in momen- 
tum. Figure IT31 fright) shows the (pt) of the momentum-balancing hadrons opposite 
to the high pt trigger as a function of centrality. The (pt) were observed to decrease 
from peripheral to central Au+Au collisions. Also shown in the Figure 11.51 (right) 
is the (pt) of the inclusive hadrons as a function of centrality. This study will be 
extended to higher pt trigger particle. The results may suggest that the moderately 
hard parton traversing a significant path length through the collision matter makes 
substantial progress toward equilibrium with the bulk. The rapid attainment of ther- 
malization via multitude of softer parton-parton interactions in the earliest collision 
stages would then not be so surprising [il3i . 

1.3.3 Particle composition in Au+Au at intermediate px 

As we have mentioned above, for Raa, the pQCD model including the parton energy 
loss, Cronin enhancement and nuclear shadowing can qualitatively fit the trend of 
data a.t 2 < Pt < Q GeV/c, however, the quantitative discrepancy between the model 
and the data is also obvious. In the intermediate pt, the mechanism for particle 
production may be different from that at high pt- 

Figure O (left) shows the vr, K, p spectra in 0%-5% and 60%-92% Au+Au 200 GeV 
collisions from ^3]. It shows that the shapes of the spectra show clear mass depen- 
dence. And in central collisions, the vr, K, p yields are close to each other at p-r > 2 
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GeV/c while it's not the case in peripheral collisions. Figure [T^ fright) shows pro- 
ton/pion (top) and anti-proton/pion (bottom) ratios for central 0-10%, mid-central 
20-30% and peripheral 60-92% in Au+Au collisions at 200 GeV HI]. It shows that 
the p{p)/7T ratios increase fast from peripheral to central collisions. In the 0-10% 
centrality bin, the proton yield is even larger than pion yield at intermediate px- 
Figure 11.71 shows the ratio Rcp for identified mesons and baryons at mid-rapidity 
calculated using centrality intervals, 0-5% vs. 40-60% of the collision cross section 
from STAR measurement ^2)- It seems that for meson, the Rcp follows a common 
trend and for baryon, the Rcp also follows a common trend, which is different from 
that for mesons. The Rcp for baryons is observed to be larger than that for mesons. 
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Figure 1.6: (left) The vr, K, p spectra in 0%-5% and 60%-92% Au+Au 200 GeV 
collisions from (right) Proton/pion (top) and anti-proton/pion (bottom) ratios 
for central 0-10%, mid-central 20-30% and peripheral 60-92% in Au+Au collisions 
at 200 GeV. Open (filled) points are for charged (neutral) pions. The data at a/s = 
53 GeV p+p collisions are also shown. The solid line is the {p + p)/ (vr+ + tt^) 
ratio measured in gluon jets jlH]. This figure is from |14j . 

These experimental results suggest that the degree of suppression depends on 
particle species(baryon/meson) at intermediate pt- The spectra of baryons (protons 
and lambdas) are less suppressed than those of mesons (pions, kaons) [HI in the 
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Pt range 2 < px < 5 GeV/c. The baryon content in the hadrons at intermediate pt 
depends strongly on the impact parameter (centrahty) of the Au+Au colhsions with 
about 40% of the hadrons being baryons in the minimum-bias colhsions and 20% 
in very peripheral collisions [HlCHl- Hydrodynamics ^21120]; parton coalescence at 
hadronization and gluon junctions have been suggested as explanations 

for the observed particle-species dependence. 
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Figure 1.7: The ratio Rcp for identified mesons and baryons at mid-rapidity calcu- 
lated using centrality intervals, 0-5% vs. 40-60% of the collision cross section. The 
bands represent the uncertainties in the model calculations of Nbm- We also show 
the charged hadron Rcp measured by STAR for ^,/s^ = 200 GeV ISlj. This figure is 
from |12j . 

In these models, recombination/coalescence models successfully reproduce Raa of 
baryons and mesons at intermediate pt, as well as showing consistency with the v2 
measurement in the same pt range. 
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Recombination model 

The concept of quark recombination was introduced to describe hadron production 
at forward rapidity in p+p collisions [23]. At forward rapidity, this mechanism allows 
a fast quark resulting from a hard parton scattering to recombine with a slow anti- 
quark, which could be one in the original sea in the incident hadron, or one incited 
by a gluon [221 • If a QGP is formed in the relativistic heavy ion collisions, then 
one might expect coalescence of the abundant thermal partons to provide another 
important hadron production mechanism, active over a wide range of rapidity and 
transverse momentum In particular, at moderate values(above the realm 

of hydrodynamics applicability), the hadron production from recombination of lower 
Pt partons from thermal bath ^Tl IIH] has been predicted to be competitive 
with the production from fragmentation of higher pt scattered partons. It has been 
suggested [2Z1 that the need for substantial recombination to explain the observed 
hadron yield and flow may be taken clS cl SI gnature of QGP formation. 
In order to explain the features of RHIC collisions, the recombination models [221 
1^ 1221 12B| make the central assumption that coalescence proceeds via constituent 
quarks, whose number in a hadron determines its production rate. The constituent 
quarks are presumed to follow a thermal (exponential) momentum spectrum and 
to carry a collective transverse velocity distributions. This picture leads to clear 
predicted effects on baryon and meson production rates, with the former depending 
on the spectrum of thermal constituent quarks and antiquarks at roughly one-third 
the baryon p^, and the latter determined by the spectrum at roughly one-half the 
meson px- Indeed, the recombination model was recently was re-introduced at RHIC 
context, precisely to explain the abnormal abundance of baryon vs meson observed at 
intermediate pt 121, J2^ ,23,. If the observed saturated elliptic flow values of hadrons 
in this momentum range result from coalescence of collectively flowing constituent 
quarks, then one expect a similarly simple baryon vs meson relationship [211 12211231: 
the baryon (meson) flow would be 3 (2) times the quark flow at roughly one-third 
(one-half) the baryon (meson) px [T^ . 
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1.3.4 Summary 

In summary, the several important results from RHIC have been introduced. The el- 
liptic flow v2 can be reproduced by hydrodynamics at low pr- At intermediate pr, v2 
from data show saturation and deviate from hydro dynamical model predictions. At 
the same time, v2 from data show baryon or meson species dependence. High px sup- 
pression can be reproduced by pQCD model and gluon saturation model. The gluon 
saturation model is also called color glass condensate model (CGC). The production 
rate dependence on baryon or meson species has been observed at intermediate Pt, 
which can be reproduced by the recombination model. 

1.4 Cronin effect 

1.4.1 Why we need d+Au run at RHIC 

In order to see the intermediate and high Pt suppression is due to the final-state effect 
or initial state effect, the measurements from d+Au collisions will provide the essential 
proof. Since the initial state in d+Au collisions is similar to that in Au+Au collisions, 
and, it's believed that the quark-gluon plasma doesn't exist in d+Au collisions, the 
results from d+Au collisions will be very important for us to judge whether the quark- 
gluon plasma exists in Au+Au collisions or not and to understand the property of the 
dense matter created in Au+Au collisions. Besides, if the identified particle spectra 
in d+Au and p+p collisions are measured, they will not only provide the reference for 
those in Au+Au collisions at 200 GeV, but also provide a chance to see the mechanism 
of the Cronin effect itself clearly at 200 GeV. Cronin effect was observed 30 years ago 
experimentally and the study of this effect was only limited to lower energy fixed 
target experiments. Before we go to the d+Au collisions, let's look back on the p+A 
collisions at lower energy fixed target experiment. 
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1.4.2 Lower energy 



The hadron spectra have been observed to depend on the target atomic weight (A) 
and the produced particle species in lower energy p+A collisions 28 . This is known 
as the "Cronin Effect", a generic term for the experimentally observed broadening 
of the transverse momentum distributions at intermediate px in p+A collisions as 
compared to those in p+p collisions j2Hl 1211 1^0] ■ The effect can be characterized as 
a dependence of the yield on the target atomic weight as A". At energies of ^/s ^ 
30 GeV, a depends on pt and is greater than unity at high pt j2H], indicating an 
enhancement of the production cross section. As shown in Figure ll.8[ the a is larger 
than 1 in the intermediate pt and shows strong particle-species dependence. The a 
for proton and antiproton are larger than those for kaon and pion. And a for kaon 
is larger than that for pion. This effect has been interpreted as partonic scatterings 
at the initial impact [211 ED] ■ Besides, the lower energy data suggest the power a 
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Figure 1.8: (left) The power alpha of A dependence from 300 GeV incident proton- 
fixed target experiment. This figure is from ^H^. (right) The Cronin ratio Rw/Be 
Pt = 4.61 GeV/c versus energy. This plot is from ,28^. 



decreases with energy, as shown in Figure 11.81 However, the energy dependence 
study of Cronin effect is limited to fixed target experiment at lower energy. What's 
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the extrapolation of Cronin effect at higher energy such as RHIC energy 200 GeV. 
At higher energies, multiple parton collisions are possible even in p+p collisions |32j. 
This combined with the hardening of the spectra with increasing beam energy would 
reduce the Cronin effect . There are several models which give different predictions 
of Cronin effect at 200 GeV. 

1.4.3 Predictions: RHIC energy 

One of the models is the initial multiple parton scattering model. In this model, the 
transverse momentum of the parton inside the proton will be broadened when the 
proton traverses the Au nucleus due to the multiple scattering between the proton and 
the nucleons inside the Au nucleus. In these models, the Cronin ratio will increases to 
a maximum value between 1 and 2 at 2.5< pt <4.5 GeV/c and then decreases with 
Pt increasing [201 ■ The Cronin effect is predicted to be larger in central d+Au colli- 
sions than in d+Au peripheral collisions [HI]. Another model is the gluon saturation 
model. At sufficiently high beam energy, gluon saturation is expected to result in a 
relative suppression of hadron yield at high pt in both p+A and A+A collisions and 
in a substantial decrease and finally in the disappearance of the Cronin effect [33J. 
Figure 11.91 shows the RdAu of charged hadron vs pr from STAR. We can see that the 
Cronin ratio increases to a maximum value around 1.5 at 3 < < 4 GeV/c and 
then decreases again jSH]- This is consistent with the initial multiple parton scatter- 
ing model [Sni- These results on inclusive hadron production from d+Au collisions 
indicate that hadron suppression at intermediate and high px in Au+Au collisions is 
due to final state interactions in a dense and dissipative medium produced during the 
collision and not due to the initial state wave function of the Au nucleus |3ti| I37j. 
Now we know that the hadron suppression at intermediate p^ in Au+Au collisions 
is due to final-state effects [211 EZl- What's the effect on particle composition at the 
same pt range in Au+Au collisions? Another question is whether there is any Cronin 
effect dependence on particle-species in d+Au collisions or not. In order to further 
understand the mechanisms responsible for the particle dependence of px spectra in 
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Figure 1.9: Rab for minimum bias and central d+Au collisions, and central Au+Au 
collisions 0. The minimum bias d+Au data are displaced 100 MeV/c to the right for 
clarity. The bands show the normalization uncertainties, which are highly correlated 
point-to-point and between the two d-l-Au distributions. This Figure is from 

heavy ion collisions, and to separate the effects of initial and final partonic rescatter- 
ings, we measured the pt distributions of vr^, K^, p and p from 200 GeV d+Au and 
p+p collisions. In this thesis, we discuss the dependence of particle production on 
Pt, collision energy, and target atomic weight. And we compare the Cronin effect of 
TT^, K^, p and p with models to address the mechanism for Cronin effect in d+Au 
collisions at ^/s^ = 200 GeV. 



Chapter 2 



The STAR 



Experiment 



2.1 The RHIC 



Accelerator 



The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Lab (BNL) is 
the first hadron accelerator and collider consisting of two independent ring. It is 
designed to operate at high collision luminosity over a wide range of beam energies 
and particle species ranging from polarized proton to heavy ion jU |321, where the 
top energy of the colliding center-of-mass energy per nucleon-nucleon pair is i/s/viv 
= 200 GeV. The RHIC facility consists of two super-conducting magnets, each with 
a circumference of 3.8 km, which focus and guide the beams. 

Figure 2.1 shows the BNL accelerator complex including the accelerators used to 
bring the gold ions up to RHIC injection energy. In the first, gold ions are acceler- 
ated to 15 MeV/nucleon in the Tandem Van de Graaff facility. Then the beam is 
transferred to the Booster Synchrotron and accelerated to 95 MeV/nucleon through 
the Tandem-to-Booster line. Then the gold ions are transferred to the Alternating 
Gradient Synchrotron (AGS) and accelerated to 10.8 GeV/nucleon. Finally they are 
injected to RHIC and accelerated to the collision energy 100 GeV/nucleon. 

RHIC's 3.8 km ring has six intersection points where its two rings of accelerating 
magnets cross, allowing the particle beams to collide. The collisions produce the 
fleeting signals that, when captured by one of RHIC's experimental detectors, pro- 
vide physicists with information about the most fundamental workings of nature. If 
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Figure 2.1: A diagram of the Brookhaven National Laboratory collider complex in- 
cluding the accelerators that bring the nuclear ions up to RHIC injection energy (10.8 
GeV/nucleon for ^^''Au). Figure is taken from |40j Slj. 
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RHIC's ring is thought of as a clock face, the four current experiments are at 6 o'clock 
(STAR), 8 o'clock (PHENIX), 10 o'clock (PHOBOS) and 2 o'clock (BRAHMS). There 
are two additional intersection points at 12 and 4 o'clock where future experiments 
may be placed [H]. 

2.2 The STAR Detector 




Figure 2.2: Perspective view of the STAR detector, with a cutaway for viewing inner 
detector systems. Figure is taken from 

The Solenoidal Tracker at RHIC (STAR) is one of the two large detector systems 
constructed at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory. STAR was constructed to investigate the behavior of strongly interact- 
ing matter at high energy density and to search for signatures of quark-gluon plasma 
(QGP) formation. Key features of the nuclear environment at RHIC are a large 
number of produced particles (up to approximately one thousand per unit pseudo- 
rapidity) and high momentum particles from hard parton-parton scattering. STAR 
can measure many observables simultaneously to study signatures of a possible QGP 
phase transition and to understand the space-time evolution of the collision process 
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Figure 2.3: Cutaway side view of the STAR detector as configured in 2001. Figure is 
taken from 



in ultra-relativistic heavy ion colhsions. The goal is to obtain a fundamental under- 
standing of the microscopic structure of these hadronic interactions at high energy 
densities. In order to accomplish this, STAR was designed primarily for measure- 
ments of hadron production over a large solid angle, featuring detector systems for 
high precision tracking, momentum analysis, and particle identification at the center 
of mass (cm.) rapidity. The large acceptance of STAR makes it particularly well 
suited for event-by-event characterizations of heavy ion collisions and for the detec- 
tion of hadron jets [IH] . 

The layout of the STAR experiment jH] is shown in Figure A cutaway side view 
of the STAR detector as configured for the RHIC 2001 run is displayed in Figure 1231 
A room temperature solenoidal magnet with a maximum magnetic field of 0.5 T 
provides a uniform magnetic field for charged particle momentum analysis. Charged 
particle tracking close to the interaction region is accomplished by a Silicon Vertex 
Tracker 03 (SVT). The Sihcon Drift Detectors gZI (SDD) installed after 2001 is also 
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for the inner tracking. The sihcon detectors cover a pseudo-rapidity range | |< 1 
with complete azimuthal symmetry (A0 = 27r). Sihcon tracking close to the interac- 
tion allows precision localization of the primary interaction vertex and identification 
of secondary vertices from weak decays of, for example, A, S, and Q. A large volume 
Time Projection Chamber EH] (TPC) for charged particle tracking and particle 
identification is located at a radial distance from 50 to 200 cm from the beam axis. 
The TPC is 4 meters long and it covers a pseudo-rapidity range | 1] \< 1.8 for tracking 
with complete azimuthal symmetry (A0 = 27r). Both the SVT and TPC contribute 
to particle identification using ionization energy loss, with an anticipated combined 
energy loss resolution (dE/dx) of 7 % (cr). The momentum resolution of the SVT and 
TPC reach a value of 5p/p = 0.02 for a majority of the tracks in the TPC. The Sp/p 
resolution improves as the number of hit points along the track increases and as the 
particle's momentum decreases, as expected 

To extend the tracking to the forward region, a radial-drift TPC (FTPC) [201 is 
installed covering 2.5 <| r] \< 4, also with complete azimuthal coverage and sym- 
metry. To extend the particle identification in STAR to larger momenta over a 
small solid angle for identified single-particle spectra at mid-rapidity, a ring imag- 
ing Cherenkov detector jM] covering | |< 0.3 and A0 = O.llvr, and a time-of-fiight 
patch (TOFp) covering — 1 < 77 < and A0 = 0.04-7? (as shown in Figure 17!^ 
was installed at STAR in 2001 ,43^. In 2003, a time-of-fiight tray (TOFr) based on 
multi-gap resistive plate chamber (MRPC) technology [^2] was installed in STAR 
detector, covering — 1 < < and A(f) = 7r/30. For the time-of-fiight system, the 
Pseudo- Vertex Position Detectors (pVPD) was installed as the start-timing detector, 
which was 5.4 m away from TPC center and covers 4.4 < \r]\ < 4.9 with the azimuthal 
coverage 19% [73] in 2003. 

The fast detectors that provide input to the trigger system are a central trigger bar- 
rel (CTB) at \ri\ < 1 and two zero-degree calorimeters (ZDC) located in the forward 
directions at ^ < 2 mrad. The CTB surrounds the outer cylinder of the TPC, and 
triggers on the fiux of charged particles in the mid-rapidity region. The ZDCs are 
used for determining the energy in neutral particles remaining in the forward di- 
rections [43J. A minimum bias trigger was obtained by selecting events with a pulse 



22 



height larger than that of one neutron in each of the forward ZDCs, which corresponds 
to 95 percent of the geometrical cross section j^ . 

2.2.1 The Time Projection Chamber 

The STAR detector I33j uses the TPC as its primary tracking device. The TPC 
records the tracks of particles, measures their momenta, and identifies the particles 
by measuring their ionization energy loss (dE/dx). Particles are identified over a mo- 
mentum range from 100 MeV/c to greater than 1 GeV/c and momenta are measured 
over a range of 100 MeV/c to 30 GeV/c jSHl. 

The STAR TPC is shown schematically in Figure 12.41 It is a volume of gas in a 
well defined uniform electric field of ~ 135 V/cm. The working gas of TPC is PIO 
gas (10% methane, 90% argon) regulated at 2 mbar above atmospheric pressure 
This gas has long been used in TPCs. Its primary attribute is a fast drift velocity 
which peaks at a low electric field. Operating on the peak of the velocity curve makes 
the drift velocity stable and insensitive to small variations in temperature and pres- 
sure [38 . The paths of primary ionizing particles passing through the gas volume are 
reconstructed with high precision from the released secondary electrons which drift 
to the readout end caps at the ends of the chamber. The drift velocity of electrons 
is 5.45 cm/yus. The uniform electric field which is required to drift the electrons is 
defined by a thin conductive Central Membrane (CM) at the center of the TPC, 
concentric field cage cylinders and the read out end caps j2H]- The readout system 
is based on Multi Wire Proportional Chambers (MWPC) with readout pads. The 
drifting electrons avalanche in the high fields at the 20 fim anode wires providing an 
amplification of 1000 to 3000. The induced charge from an avalanche is shared over 
several adjacent pads, so the original track position can be reconstructed to a small 
fraction of a pad width. There are a total of 136,608 pads in the readout system [3H| , 
which give x-y coordinate information. The z position information is provided by 512 
time buckets. 

At the Data Acquisition (DAQ) stage, raw events containing millions of ADC values 
and TDC values were recorded. Raw data were then reconstructed into hits, tracks. 
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Figure 2.4: The STAR TPC surrounds a beam-beam interaction region at RHIC. The 
coUisions take place near the center of the TPC. 

vertices, and the colhsion vertex through the reconstruction chain of TPC by 
Kalman method. The colhsion vertex are called the primary vertex. The tracks are 
called the global tracks. If the 3-dimensional distance of closest approach (DCA/dca) 
of the global track to the primary vertex is less than 3 cm, this track will be chosen 
for a re-fit by forcing a new track helix ending at the primary vertex. These newly 
reconstructed helices are called primary tracks jHI]. As expected, the vertex resolu- 
tion decreases as the square root of the number of tracks used in the calculation. The 
vertex resolution is 350 fim when there are more than 1,000 tracks ^38;. Figure 
shows the beam's eye view of a central Au+Au collision event in the STAR TPC. 
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Figure 2.5: Beam's eye view of a central Au+Au collision event in the STAR Time 
Projection Chamber. This event was drawn by the STAR online display. Figure is 
taken from [IH] . 

Particle Identification (PID) of TPC by dE/dx 

Energy lost in the TPC gas is a valuable tool for identifying particle species. It works 
especially well for low momentum particles but as the particle energy rises, the energy 
loss becomes less mass- dependent and it is hard to separate particles with velocities 
V > 0.7c [33- For a particle with charge z (in units of e) and speed (3 = v/c passing 
through a medium with density p, the mean energy loss it suffers can be described 
by the Bethe-Bloch formula 

(^) = 27riVor,meC P^[ln 2/3 ] (2.1) 

where Nq is Avogadro's number, rrie is the electron mass, (= /m^) is the classical 
electron radius, c is the speed of light, Z is the atomic number of the absorber, 
A is the atomic weight of the absorber, 7=1/ a/1 — / is the mean excitation 
energy, and Em (= 2meC^(3'^ / (1 — is the maximum transferable energy in a single 
collision [HHll^- From the above equation, we can see that different charged particles 
(electron, muon, pion, kaon, proton or deuteron) with the same momentum p passing 
through the TPC gas can result in different energy loss. Figure ESI shows the energy 
loss for particles in the TPC as a function of the particle momentum, which includes 
both primary and secondary particles. We can see that charged pions and kaons 
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can be identified up to about transverse momentum 0.75 GeV/c and protons and 
anti-protons can be identified to 1.1 GeV/c. 
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Figure 2.6: The energy loss distribution for primary and secondary particles in tfie 
STAR TPC as a function of tlie pt of tlie primary particle. This figure is taken 
from 



In order to quantitatively describe the particle identification, we define the variable 
NcTTT (in the case of charged pion identification) as 

,0.55 ciE 



dE 

dx meas. 



dx ■ 



(2.2) 



in which is the number of hits for a track in the TPC, 4^ is the measured 

' dx meas. 

energy loss of a track and {^)n is the mean energy loss for charged pions. In order 
to identify charged kaons, protons and anti-protons, we can have similar definition of 
Na-K and N^p. Thus we can cut on the variables Na-^, N„k and N„p to select different 
particle species [8T]|. 

A specific part of the particle identification is the topological identification of neutral 
particles, such as the Kg and A. These neutral particles can be reconstructed by iden- 
tifying the secondary vertex, commonly called VO vertex, of their charged daughter 



decay modes. Kg 



and A 



pn 



8J. 
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2.2.2 The time-of-flight tray based on MRPC technology 




Figure 2.7: Tray structure. Figure is taken from |74j . 

In 2003, the time-of-flight tray (TOFr) based on muhi-gap resistive plate chamber 
(MRPC) technology jSHl was installed in STAR detector. It extends particle identi- 
fication up to ~ 3 GeV/c for p and p. This tray was installed on the Au beam 
outgoing direction. MRPC technology was first developed by the CERN ALICE 
group [SEj to provide a cost-effective solution for large-area time-of-flight coverage. 
For full time-of-flight coverage at STAR, there will be 120 trays, with 60 on east 
side and 60 on west side. For each tray, there will be 33 MRPCs. For each MRPC, 
there are 6 read-out channels. Figure 12.71 shows the tray which indicates the position 
of each MRPC module. The MRPCs are tilted differently so that each MRPC is 
most projective to the average primary vertex location at Z=0. In 2003 d+Au and 
p+p run, only 28 MRPCs were installed in the tray and 12 out of 28 were instru- 
mented with the electronics, representing 0.3% of TPC coverage. If we number the 
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33 MRPCs in the tray from 1 to 33, with 1 close to TPC center and 33 far from 
TPC center, the numbers of 12 modules instrumented with the electronics in 2003 
are 3,4,5,7,9,10,11,12,13,14,26 and 32. 



The introduction of MRPC 
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Figure 2.8: Two side views of MRPC. The upper (lower) is for long (short) side view. 
The two plots are not at the same scale. Figure is taken from [71]. 



Resistive Plate Chambers (RPCs) were developed in 1980s [HZ], and were originally 
operated in streamer mode. This operation mode allows us to get high detection 
efficiency (>95%) and time resolution ( 1 ns), with low fluxes of incident particles. At 
higher fluxes (>200 Hz/cm?)^ RPCs begin to lose their efficiency. A way to overcome 
this problem is to operate RPCs in avalanche mode. The Multi-gap Resistive Plate 
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Figure 2.9: The shape of the 6 read-out strips for each MRPC. 

Chamber (MRPC) was developed less than 10 years ago [5H]. It consists of a stack of 
resistive plates, spaced one from the other with equal sized spacers creating a series 
of gas gaps. Electrodes are connected to the outer surfaces of the stack of resistive 
plates while all the internal plates are left electrically floating. Initially the voltage 
on these internal plates is given by electrostatics, but they are kept at the correct 
voltage due to the flow of electrons and ions created in the avalanches. Figure 12.81 
shows the structure of MRPC detector. For each MRPC, there are 6 read-out strips. 
Figure 12.91 shows the shape of the read-out strip. The detailed production process 
can be found at Appendix B. 

MRPC, as a new kind of detector for time of fight system, operated in avalanche mode 
with a non flammable gas mixture of 90% F134A, 5% isobutane, 5% SF6, can fulfill 
all these requirements: high efficiency (>95%), excellent intrinsic time resolution 
(<100 ps) inniEllZnilZIlIZS, high rate capability ( 500 Hz/cm'^), high modularity 
and simplicity for construction, good uniformity of response, high granularity/ low 
occupancy, and large acceptance. 

Simulation: the work principle of this chamber 

A detailed description of the model used in the simulation was reported in these pa- 
pers 1211 inni EH E21 , here just the main items will be repeated. The program starts 
from considering an ionizing particle which crosses the gas gaps and generates a cer- 
tain number of clusters of ion-electron pairs. The electrons contained in the clusters 
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drift towards the anode and, if the electric field is sufficiently high, give rise to the 
avalanche processes. 

The primary cluster numbers and the avalanche growth are assumed to follow, re- 
spectively, simple Poisson statistics and the usual exponential law. Avalanche gain 
fluctuations have been taken into account using a Polya distribution After the 
simulation of the drifting avalanches, the program computes, by means of Ramo [Ul] 
theorem, the charge qind induced on the external pick-up electrodes (strips or pads) 
by the avalanche motion. Under certain approximations, this is given by the formula 



where qe is the electron charge, t] 1st effective Townsend coefficient rj = a — (3, a 
is the Townsend coefficient, (3 is the attachment coefficient, Xj the jth cluster initial 
distances from the anode, d the gap width, rij the number of initial electrons in the 
considered jth cluster, M the avalanche gain fluctuations factor, and AVw/d = is 
the normahzed weighting field. In addition to qind, the current iind{t) induced on the 
same electrodes by the drifting charge qd{t) may be computed as 



, where Vd is the electron drift velocity. The computation of imd allows us to repro- 
duce the whole information coming out from MRPC, such as time distribution. 
Charge Spectrum Simulation: The almost Gaussian charge distribution obtained 
with the MRPC is a key ingredient to its performance. If the avalanches grew follow- 
ing Townsend's formula the charge distribution would be exponential in shape. Thus 
the space charge effects must be considered in the simulation. 

The input parameters for the simulation program are: the Townsend coefficient a, the 
attachment coefficient /5, the average distance between clusters A and the probability 
distribution of the number of electrons per cluster. These pieces of information can 
be obtained, for a given gas mixture and given conditions (pressure and temperature) 
and electric field, by the programs HEED and MAGBOLTZ [Ml I67j . In addition, 
a maximum number of electrons in an avalanche (cutoff value) is specified. 




(2.3) 



iiUt) = AV^^qdit)Me 



(2.4) 
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In a given gap, we generate a number of clusters with distances exponentially dis- 
tributed with average distance A. For each cluster, we then generate a certain num- 
ber of electrons, according to the distribution obtained by the program HEED. Each 
electron from the primary cluster will give rise to a number of electrons, generated 
according to an exponential probability law. 
For each cluster, the avalanche growth is stopped when the total charge reaches 




Figure 2.10: Simulated 1st effective Townsend coefficient curve and normalized charge 
distribution for a 6 and 10 gap MRPC 

a certain cutoff value, as originally suggested in ref. [HH] to take into account space 
charge effects in the avalanche development. This cutoff value has been set to be 
1.6 X 10^ electrons. 

In Figure 12.101 we show the results of simulations. Figure I2.10f a) is the simulated 
curve of the 1st effective Townsend coefficient rj versus the electric field, which is 
generated by Magboltz. The curve shows that the correlation between rj and the 
electric field is almost linear when MRPC is operated at high electric field for the gas 
mixture. Figure O.lOf b) is the charge spectrum for a 6 gap chamber and (c) (d) for a 
10 gap chamber compared to experimental data jnHEH]; and the number under each 
plot shows the electric field E in the gas gap for MRPC. In both cases the gap size 
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is 220 fim. The gas mixture was 90% F134A, 5% isobutane and 5% SF6 in normal 
conditions of pressure and temperature. The value of A used was 0.1 mm, derived 
from HEED program. 

The charge distribution has an almost Gaussian form, especially for the 10 gap MRPC 
The left side of the distribution (very few events at values near zero) is due to the 
fact that the MRPC operates at high gain r/ x d ~ 30. This means that avalanches 
starting in the middle of the gap width, which only avalanche over half the distance, 
give a detectable signal. The charge distribution is the superposition of several prob- 
ability distributions which, according to the central limit theorem, will tend to a 
Gaussian form. The right side of the charge distribution (the fact that the tails are 
not very long) indicates that indeed the space charge effects stop the development of 
the avalanche. 

Time Distribution Simulation: We then proceed to simulate the time distribu- 
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Figure 2.11: Simulated results of a 6 gap MRPC. 

tion of these same chambers. The electron drift velocity can be obtained from HEED. 
When the total induced charge signal is over threshold, the time is recorded. In this 
paper, the threshold is 13 fc for the 6 gap MRPC and 26 fc for the 10 gap MRPC. 
Fig. 2 is the simulated results for a 6 gap chamber. Figure 12.111 (a) is the simulated 
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curve of the electron drift velocity versus the electric field, which is generated by Mag- 
boltz. Figure 12.111 (b) is the time distribution of a 6 gap MRPC The intrinsic time 
resolution is only 19 ps or so. If we consider other contributions, such as front-end 
electronics 30 ps, TDC resolution 25 ps, fanout start signal 10 ps, beam size (1cm) 15 
ps, we can get the MRPC resolution is ^20^ + 30^ + 25^ + 10^ + 15^ = 47 ps. This 
value is similar to the experimental result [021 EHI- For a 10 gap MRPC, the intrinsic 
time resolution is about 15 ps. 

From the simulation, we can get the bottom line of MRPC time resolution ~ 20 ps. 
And we need to keep control of all these contributions to ensure best time resolution. 

MRPC for this tray installed in 2003 

In 2003, for the MRPCs in the TOFr, the inner glass thickness is 0.54 mm, the outer 
glass is 1.1 mm. The gas gap is 0.22 mm. Both the volume resistivity {10^'^~^^ohm.cm) 
of the glass plates and the surface resistivity(2M ohm per square) of carbon layer at 
room temperature are presented in [7S]. It is found the volume resistivity of the 
plate decreases with the temperature increasing. And the radiation will decrease the 
volume resistivity of the plate [73]. In order not to pollute the working gas of TPC, 
SF6 is not used as part of the working gas of TOFr. The working gas of MRPC-TOFr 
at STAR is 95% freon and 5% iso-butane at normal atmospheric pressure. The high 
voltage applied to the electrodes is 14.0 kV. 



Chapter 



3 



Analysis 



Methods 



3.1 Trigger 



The detector used for these studies was the Solenoidal Tracker at RHIC (STAR). The 
main tracking device is the Time Projection Chamber (TPC) which provides mo- 
mentum information and particle identification for charged particles up to ~ 1-1 
GeV/c by measuring their ionization energy loss (dE/dx) [38 . Detailed descriptions 
of the TPC and d+Au run conditions have been presented in Ref. j36| I38j. A pro- 
totype time-of-fiight detector (TOFr) based on multi-gap resistive plate chambers 
(MRPC) 39 was installed in STAR for the d+Au and p+p runs. It extends particle 
identification up to ~ 3 GeV/c for p and p. 

TOFr covers 7r/30 in azimuth and — 1 < r/ < in pseudorapidity at a radius of 
~ 220 cm. It contains 28 MRPC modules which were partially instrumented during 
the 2003 run. Since the acceptance of TOFr is small, a special trigger selected events 
with a valid pVPD coincidence and at least one TOFr hit. A total of 1.89 million 
and 1.08 million events were used for the analysis from TOFr triggered d+Au and 
non-singly diffractive (NSD) p+p collisions, representing an integrated luminosity of 
about 40 /ib~^ and 30 nb~^, respectively. Minimum-bias d+Au and p+p collisions 
that did not require pVPD and TOFr hits were also used to study the trigger bias 
and enhancement, and the TOFr efficiency and acceptance. The d+Au minimum-bias 
trigger required an equivalent energy deposition of about 15 GeV in the Zero Degree 
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Calorimeter in the Au beam direction The trigger efficiency was determined 

to be 95 ± 3%. Minimum-bias p+p events were triggered by the coincidence of two 
beam-beam counters (BBC) covering 3.3 < \r]\ < 5.0 [8^. The NSD cross section was 
measured to be 30.0 ± 3.5 mb by a van der Meer scan and PYTHIA [7^1 simulation 
of the BBC acceptance 0. 

3.1.1 Centrality tagging 

Centrality tagging of d+Au collisions was based on the charged particle multiplicity 
in —3.8 < f] < —2.8, measured by the Forward Time Projection Chamber in the 
Au beam direction 175] . The TOFr triggered d+Au events were divided into 
three centralities: most central 20%, 20 — 40% and 40— ~ 100% of the hadronic 
cross section. The average number of binary collisions (Nhin) for each centrality class 
and for the combined minimum-bias event sample is derived from Glauber model 
calculations and listed in Table 13.11 

Table 13.11 also lists the uncorrected FTPC east reference multiplicity ranges for 
centrality definitions. 



Centrality Bin 


Uncorr. FTPCRefMult Range 


Uncorr. N^harge 




M.B. 




10.2 


7.5 ±0.4 


0%-20% 


FTPCRefMult > 17 


17.58 


15.0 ± 1.1 


20%-40% 


10 < FTPCRefMuh < 17 


12.55 


10.2 ± 1.0 


40%-100% 


< FTPCRefMuh < 10 


6.17 


4.0 ±0.3 



Table 3.1: Centrality definitions for different uncorrected FTPC east reference multi- 
plicity ranges. Uncorrected Ncharge stands for the average value of uncorrected refer- 
ence multiplicity in certain centrality bin. The fourth column represents the number 
of binary coUisons {Nun) calculated from Glauber model. 



3.1.2 Trigger bias study 

Since we set up a special trigger which selected events with a valid pVPD coincidence 
and at least one TOFr hit, the study of pr dependence of trigger bias is neces- 
sary. Figure Em shows there is negligible trigger bias on pt dependence at > 0.3 
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GeV/c from simulation. In this figure, pVPD means that pVPD is required to fire in 
minimum-bias colhsions. TOF means that TOFr is required to fire in minimum-bias 
coUisions, and pVPD & TOF means that pVPD and TOFr are required to fire in 
minimum-bias coUisions. From this figure, if we required pVPD and TOFr to fire, 
we can see the ratio is flat with px when pt is larger than 0.3 GeV/c by comparison 
through the px distribution in minimum-bias collisions. That means the trigger bias 
for Pt distribution is negligible at p-p > 0.3 GeV/c. 



tracks in TOF acceptance 




0.5 1 1.5 2 2J 3 3.5 4 4^ 5 

Pr fGeVfa) 

Figure 3.1: The px dependence plot of the trigger bias. 



Minimum-bias d+Au and p+p collisions are used to study the trigger bias and 
enhancement. Figure 021 shows the trigger bias and enhancement in d+Au minimum- 
bias collisions and three centrality bins. In this figure, TOFr means that TOFr is 
required to fire in minimum-bias events. pVPD means that TOFr and pVPD are 
required to fire in minimum-bias events. Minbias means the minimum-bias triggered 
events. For enhancement study, TOFr/pVPD is the ratio of the number of events in 
which TOFr is required to fire over the number of events in which TOFr and pVPD 
are required to fire, and Minbias/pVPD is the ratio of the number of minimum- bias 
triggered events over the number of events in which TOFr and pVPD are required to 
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Figure 3.2: The enhancement factor and (Nch) bias in minimum-bias and centrahty 
selected d+Au coUisions. 

fire. The enhancement factor for TOFr is (Minbias/pVPD)/(TOFr/pVPD). For ex- 
ample, in minimum-bias collision, Minbias/pVPD is equal to 28.7, while TOFr/pVPD 
is 2.87, so in minimum-bias collisions, the enhancement of TOFr trigger is 10. For 
(Nch) bias study, TOFr/pVPD is the ratio of {Net) in the events where TOFr is re- 
quired to fire over the (Net) in the events where TOFr and pVPD are required to 
fire. Since in our triggered events, TOFr and pVPD are required to fire, TOFr/pVPD 
is our (Nch) bias factor. The curves in this figure show the charged particle multi- 
plicity at mid-rapidity in TOFr events and in TOFr and pVPD events individually. 
Table 13.21 lists the enhancement factor and trigger bias in minimum-bias, centrality 
selected d-l-Au collisions and minimum-bias p+p collisions. 



3.2 Track selection and calibration 

The TPC and TOFr are two independent systems. In the analysis, hits from par- 
ticles traversing the TPC were reconstructed as tracks with well defined geometry, 
momentum, and dE/dx The particle trajectory was then extended outward to 
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Centrality Bin 


TOFr triggered events 


enhancement factor 


{Nch) bias 


0%-100% 


1.80 M 


10.0 


1.02 


0%-20% 


0.523 M 


5.75 


1.04 


2{fA-m 


0.500 M 


8.03 


1.03 


40%-100% 


0.479 M 


15.8 


0.965 


p+p 


0.995 M 


37.4 


1.19 



Table 3.2: Trigger bias study. The {Nch) bias and enhancement factor in minimum- 
bias, centrality selected d+Au collisions and minimum-bias p-|-p collisions. 

the TOFr detector plane. The pad with the largest signal within one pad distance 
to the projected point was associated with the track for further time-of-flight and 
velocity calculations. 

3.2.1 Calibration 

pVPD calibration 

For TOFr, we use pVPD as our start-timing detector. In d-l-Au and p-|-p collisions, 
at least one east pVPD and one west pVPD were required to fire. In d-l-Au collisions, 
to cahbrate east pVPD, we required 3 east pVPD to fire; to calibrate west pVPD, we 
required 3 west pVPD to fire. In p-|-p coUisions, to cahbrate east pVPD, we required 
2 east pVPD to fire; to calibrate west pVPD, we required 2 west pVPD to fire. Let's 
take the east pVPD calibration in d-l-Au coUisions as an example. The label for 3 
pVPD are pVPDl, pVPD2, pVPD3, the adc and tdc value for pVPDl are al, tl, and 
the slewing correction function is /I; the adc and tdc value for pVPD2 are a2, t2, 
and the slewing correction function is /2; the adc and tdc value for pVPD3 arc a3, 
t3, and the slewing correction function is /3. We use tl — ((t2 — /2) + (t3 — /3))/2 vs 
al to get the slewing correction for pVPDl; use t2 — ((t3 — /3) + (tl — /l))/2 vs a2 
to get the slewing correction for pVPD2; use t3 — ((tl — /I) + (t2 — /2))/2 vs a3 to 
get the slewing correction for pVPD3. At the beginning, /I = /2 = /3 = 0, wc got 3 
curves of tl - ((t2 - /2) + (t3 - /3))/2 vs al, t2 - ((t3 - /3) + {tl - fl))/2 vs a2 and 
t3— ((tl — /l) + (t2— /2))/2 vs a3. The 3 curves corresponded to the 3 slewing functions 
/(al),/(a2),/(a3); For the second step, fl = /(al),/2 = /(a2),/3 - /(a3), also 
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plot tl - ((t2 - /2) + {tS - /3))/2 vs al, t2 - ((t3 - /3) + (tl - /l))/2 vs a2 and 
t3 — ((tl — /I) + (t2 — /2))/2 vs a3. And we got the new three slewing curves 
/'(al),/'(«2),/'(a3). For the third step, /I = /'(al),/2 = /'(a2),/3 = /'(a3), also 
plot tl - ((t2 - /2) + (t3 - /3))/2 vs al, t2 - ((t3 - /3) + (tl - /l))/2 vs a2 and 
t3 — ((tl — /I) + (t2 — /2))/2 vs a3. And we got another new three slewing curves 
/"(al), /"(a2), /"(a3). And so on and so forth till the resolution of tl-/l-((t2-/2) + 
(t3 - /3))/2, t2 - /2 - ((t3 - /3) + (tl - /l))/2 and t3 - /3 - ((tl - /I) + (t2 - /2))/2 
converged. The looping method is to subtract the correlation of different pVPD 
tubes in the same direction. The function for the slewing correction we use is y = 
par[0] + par[l]/ y/x + par[2]/x + par[3] x x. In Figure ESI the left plot shows the 
pVPD2 slewing plot and the right plot shows that the timing is independent on the 
ADC value after the slewing correction. 
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Figure 3.3: pVPD slewing correction. 

After the slewing correction, we got the corrected timing of east pVPD and west 
pVPD. For each side, the timing difference should be shifted to zero. That's to say 
the mean value in the distribution of tl — /I — (t2 — /2) and tl — /I — (t3 — /3) 
were shifted to zero. Also we need to correct for the effect caused by the different 
numbers of fired pVPD in different events. What we did was shifting the mean value 
of the distribution of {Ylite)/Ne — {^tw)/Nw — 2. x Vz/c to zero, where the ^te, 
^ tw means the sum of the corrected timing of east fired pVPD and west fired pVPD 
respectively, iVe, Nw means the number of east fired pVPD and west fired pVPD, Vz 
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is the z value of primary vertex of the event, and c is the hght velocity. 



TOFr calibration 

After the slewing correction for pVPD, we use this variable as our start timing: 



start 
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Figure 3.4: dE/dx vs p plot from d+Au collisions. The line represents that dE/dx = 
0.028 X 10^^ GeV/cm in this momentum range 0.3 < p < 0.6 GeV/c. 

The difference between TOFr timing Tt^jr and start timing Tgtart is our time of 
flight tof = Ttofr — T start- To Calibrate the to/, the pure pion sample was chosen by 
selecting the particle energy loss dE/dx in TPC at dE/dx < 0.028 x 10~^ GeV/cm in 
the momentum range 0.3 < p < 0.6 GeV/c. Figure EiU shows dE/dx vs p plot from 
d+Au collisions. Firstly the so called Tq correction was done due to the different cable 
lengths for different read-out channels, which was done by shifting the mean value of 
the distribution of tof — T^r to zero channel by channel, where T^r is the calculation 
timing assuming the particle was pion particle. Secondly, the slewing correction due 
to correlation between timing and signal amplitude of the electronics was done by 
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Figure 3.5: The slewing correction. 

getting the curve of tof — vs adc for each channel, where the tof was the time 
of flight after the Tq correction and adc was the ADC value of TOFr. The slewing 
curve is like the plot shown in Figure 13.51 The function of the slewing correction is 
y = par[0] + par[l]/^/x + par[2]/x + par[3]/ y/x/x + par[A\/x/x. 

The z position correction was also done since the different hit positions on the read- 
out strip will generate different transmission timing. This was done by getting the 
function of to/" — versus Ziocai, where the tof" is the time of flight after the Tq and 
slewing correction, and Ziocai is the the hit local z position of the TOFr. The function 
for the z position correction is y = Yli=o x 2;*)- The z position correction for 

all the channels is shown in Figure ESI After the z position was done, the calibration 
for TOFr was finished. The overall resolution of TOFr was 120 ps and 160 ps in d+Au 
and p+p collisions respectively, where the effective timing resolution of the pVPDs 
was 85 ps and 140 ps, respectively. Figure 13.71 shows the overall resolution of TOFr 
in d+Au and p+p collisions. 
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Figure 3.6: The z position correction. 

3.3 Raw yield 

Prom the timing information t from TOFr after the cahbration and the pathlength L 
from TPC, the velocity /5 of the particle can be easily got hj j3 = L/t/c. Figure ESI 
shows 1//3 from TOFr measurement as a function of momentum (p) calculated from 
TPC tracking in TOFr triggered d+Au collisions. The raw yields of vr^, K"^, p and p 
are obtained from Gaussian fits to the distributions in = p^(l//5^ — 1) in each px 
bin. 

3.3.1 TT raw yield extraction 

For TT^, the rapidity range is —0.5 < y.,, < 0.. After |A''o-7r| < 2 was required, the 
mass squared = p^(l//5^ — 1) distributions in different px bin in d+Au minimum- 
bias collisions are shown is Figure ITTHl and Figure 13'. 171 At px < 0.8 GeV/c, the 
single Gaussian function was used to fit the distribution of to get the raw yield. 
At the same time, the counting result by counting the track number at the range 
—0.1 < < 0.1 {GeV/c^)"^ was also used to compare with the raw yield from the 
fitting method. The difference between them was found in one sigma range. The raw 
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Figure 3.7: The overall timing resolution after the calibration. 

yield we quote is from the fitting method. At pr > 0.8 GeV/c, the double Gaussian 
function was used to extract the raw yield. The raw signals in each Pt bin are shown 
in Table 13.31 and Table 13.41 Also shown in the tables are those in centrality selected 
d+Au collisions and minimum-bias p+p collisions. 



3.3.2 K raw yield extraction 

For K"^, the rapidity range is —0.5 < yx < 0. After {N^jkI < 2 was required, the 
mass squared = ( 1 — 1 ) distributions in different px bin in d+ Au minimum- 
bias collisions are shown is Figure EHH and Figure 13'. 191 At px < 0.8 GeV/c, the 
single Gaussian function was used to fit the distribution of to get the raw yield. 
At the same time, the counting result by counting the track number at the range 
0.16 < < 0.36 {GeV/c^)"^ was also used to compare with the raw yield from the 
fitting method. The difference between them was found in one sigma range. The raw 
yield we quote is from the fitting method. At pr > 0.8 GeV/c, the double Gaussian 
function was used to extract the raw yield. The raw signals in each Pt bin are shown 
in Table 13.51 and Table 13.61 Also shown in the tables are those in centrality selected 
d+Au collisions and minimum-bias p+p collisions. 
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p (GeV/c) 

Figure 3.8: 1/(3 vs. momentum for vr^, K"^, and p{p) from 200 GeV d+Au collisions. 
Separations between pions and kaons, kaons and protons are achieved up to ^ 1.6 
and 3.0 GeV/c, respectively. The insert shows m? = — 1) for 1.2 < pt < 1.4 

GeV/c. Clear separation of vr, K and p is seen. 

3.3.3 p and p raw yield extraction 

For p, the rapidity range is —0.5 < yp < 0. After |A^o-p| < 2 was required, the 
mass squared = p^(l//3^ — 1) distributions in different pt bin in d+Au minimum- 
bias colhsions are shown is Figure I3'.21l At px < 1.6 GeV/c, the single Gaussian 
function was used to fit the distribution of to get the raw yield. At the same time, 
the counting result by counting the track number at the range 0.64 < < 1.44 
{GeV/c^y was also used to compare with the raw yield from the fitting method. The 
difference between them was found in one sigma range. The raw yield we quote is 
from the fitting method. At pt > 1.6 GeV/c, the double Gaussian function was used 
to extract the raw yield. The raw signals in each Pt bin are shown in Table 13.81 
For the p, the raw yield extraction method is the same as p except that at pt < 1-Q 
GeV/c, we use the method Np = Npdca<i.cm x {Npdca<3.cm/Npdca<i.cm) to reject the 
background, where Np and Np are the number of the p and p tracks individually, and 
Npdca<i.cm/Npdca<3.cm IS the ratio of p tracks at dca < 1.0 cm over those at dca < 3.0 
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Figure 3.9: the ratio of p at dca < 1.0 cm over p at (ica < 3.0 cm. 



cm. In Figure 15.201 the first 10 px bins are for dca < 1.0 cm, the last 4 pt bins are 
for dca < 3.0 cm. Figure ITTIl shows the ratio of p tracks at dca < 1.0 cm over those 
at dca < 3.0 cm. After this correction of Np = Npdca<i.cm x {Npdca<3.cm/Npdca<i.cm), 
the p raw signals in each Pt bin are shown in Table 13.71 



3.4 Efficiency and acceptance correction 

Acceptance and efficiency were studied by Monte Carlo simulations and by matching 
TPC track and TOFr hits in real data. TPC tracking efficiency was studied by Monte 
Carlo simulations. The simulated vr^, K^, p and p are generated using a flat pr and 
a fiat y distribution and pass through GSTAR [S^ (the framework software package 
to run the STAR detector simulation using GEANT [H Ei) and TRS (the TPC 
Response Simulator [H2])- The simulated vr^, K^, p and p are then combined with 
a real raw event and we call this combined event a simulated event. This simulated 
event is then passed through the standard STAR reconstruction chain and we call 
this event after reconstruction a reconstructed event. The reconstructed information 
of those particles in the reconstructed event is then associated with the Monte-Carlo 
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Figure 3.10: TPC reconstruction efficiency of vr^, K^, p and p as a function of pr- 
The left plot for charged plus particle and the right for charged minus particle. 

information in the simulated event. And then we get the total number of simulated 
TT^, K^, p and p from simulated events in a certain transverse momentum bin. Also 
we can get the total number of associated tracks in the reconstructed events in this 
transverse momentum bin jHI]. In the end, take the ratio of the number of associated 
TT^, K^, p and p over the number of simulated vr^, K^, p and p and this ratio is 
the TPC reconstruction efficiency for a certain transverse momentum bin in the mid- 
rapidity range. Figure I3.1UI shows the TPC reconstruction efficiency of tt^, K^, p 
and p as a function of pt- 

The Matching Efficiency from TPC to TOFr were studied in real data, and the 
formula are 

j^jj TofrMatchedTracks/dAuTOFr Events (3 2) 

{MinBiasTracks / MinBiasEvents)pv PD x factorl x factor2 
where the TofrMatchedTracks/dAuTOFr Events is the number of TOFr matched 
tracks per dAuTOFr trigger event, [MinBiasTracks / MinBiasEvents)pv pd is the 
number of minimum-bias tracks per minimum-bias event by requiring the pVPD 
to &e, factorl is the enhancement factor of dAuTOFr trigger, and factor2 is the 
other factors such as the TOFr trip factor. The E f f Match includes the detector 
response efficiency. Figure 13.111 shows the matching efficiency of different particle 
species including the detector response versus px- The detector response efficiency, 
including the material absorption and scattering effect between TPC and TOFr, as 
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Figure 3.11: Matching efficiency from TOFr to TPC of vr^, K^, p and p as a function 
of pt, including detector response. The left plot for charged plus particle and the 
right for charged minus particle. 



a function of pt is shown in Figure 13.121 which is around 90% at pt > 0.3 GeV/c. 
After the material absorption and scattering effect correction, the detector response 
efficiency is around 95%. 



3.5 Background correction 

Weak-decay feeddown (e.g. vr+vr") to pions is ~ 12% at low px and ~ 5% at 

high Pt, and was corrected for using PYTHIA [TUj and HIJING |73 simulations, as 
shown in Figure 13.131 For tt spectra, the n misidentification was also corrected for, 
which is also shown in Figure 13. 131 

Inclusive p and p production is presented without hyperon feeddown correction. 
p and p from hyperon decays have the same detection efficiency as primary p and 
p [IBs contribute about 20% to the inclusive p and p yield, as estimated from the 
simulation. However, for p, there is still some scattering contribution which comes 
from the beam pipe interaction after the cut of dca < 1.0 cm. Figure 13. 141 shows the 
contribution of scattering effect for proton when we cut dca < 1.0 cm. The correction 
is done at pr < 1-1 GeV/c and negligible at higher pt- 
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Figure 3.12: The TOFr response efficiency as a function of pr- 

3.6 Energy loss correction 

The energy loss effect due to the interaction with the detector material was also 
corrected for. This was studied by simulation. Figure IH. 151 shows the momentum and 
transverse momentum correction for energy loss effect. At px >0.35 GeV/c, for vr, the 
energy loss effect is negligible while for kaon and proton, the energy loss correction 
is non-negligible at lower pj- and negligible at higher p^- The correction was done by 
shifting the position of px in the px spectra. 

3.7 Normalization 

The efficiency including vertex efficiency and trigger efficiency is 91% in d+Au minimum- 
bias collisions and 85% in p+p and 40-100% d+Au collisions. In 0%-20% and 20%- 
40% d+Au collisions, the efficiency is 100%. Since the statistic of p+p minimum-bias 
events in run 3 is not good enough for us to get very precise enhancement factor and 
Nch bias factor. We compare the vr spectra in the ffist 5 pt bin with those from the 
paper and get the additional normalization factor for p+p collisions. 
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Figure 3.13: tt background contribution as a function oipT- The circled symbols repre- 
sent the total TT background contribution including feed-down and fi misidentification. 
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contributions individually. 
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Figure 3.14: The p scattering effect contribution when we cut dca < 1.0 cm. 
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Figure 3.15: (left) p energy loss correction of different particle species as a function 
of p. Prec is the reconstructed momentum before the energy loss correction, pmc is 
the momentum after energy loss correction from simulation, 6 is the angle between 
the reconstructed momentum and beam line, (right) Pt energy loss correction of 
different particle species as a function of pt- prif^Gc) is the reconstructed transverse 
momentum before the energy loss correction, pt{MC) is the transverse momentum 
after energy loss correction from simulation. 



PT (GeV/c) 


d+Au Trigger 


0%-20% 


20%-40% 


40%-100% 


P+P 


0.3-0.4 


2.929e -1- 04 ± 171.4 


9219 ±96.21 


8735 ± 93.6 


8604 ± 92.84 


1.806e + 04 ± 134.4 


0.4-0.5 


2.185e -1- 04 ± 147.8 


6894 ± 83.03 


6657 ±81.59 


6325 ± 79.53 


1.274e + 04 ±114.3 


0.5-0.6 


1.592e -1- 04 ± 126.2 


5162 ± 71.85 


4901 ± 70.3 


4534 ± 67.34 


9180 ±95.81 


0.6-0.7 


1.166e-|-04± 108 


3832 ±62.19 


3556 ± 59.64 


3311 ± 57.54 


6531 ± 80.82 


0.7-0.8 


8556 ±92.5 


2909 ± 53.93 


2628 ± 51.26 


2368 ± 48.67 


4447 ± 66.74 


0.8-0.9 


6198 ±78.86 


2099 ± 45.85 


1936 ± 44.33 


1693 ±41.17 


2973 ± 54.57 


0.9-1 


4520 ±67.25 


1487 ±38.57 


1361 ±36.9 


1276 ± 35.74 


2132 ±46.31 


1-1.1 


3312 ±57.61 


1147 ±33.9 


1033 ±32.17 


845.9 ± 29.15 


1386 ± 37.71 


1.1-1.2 


2406 ± 49.35 


788.6 ±28.19 


752.5 ±27.58 


652.7 ±25.7 


959.2 ± 31.93 


1.2-1.4 


3227 ± 58.17 


1132 ±34.28 


934.4 ±30.98 


831.5 ± 29.82 


1183 ±40.11 


1.4-1.6 


1756 ±45.2 


573.7 ± 26.2 


543.7 ± 24 


412.8 ± 21.83 


625.5 ± 30.16 


1.6-1.8 


1046 ± 39 


337.9 ± 20.42 


309.8 ±18.62 


234 ± 16.58 


364.5 ± 32.64 



Table 3.3: 7r+ raw signal table in minimum-bias, centrality selected d+Au collisions 
and minimum-bias p+p collisions. 
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Figure 3.16: 7t~^ raw yields versus mass squared distribution. The histograms are our 
data. The curves are Gaussian fits. 
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Figure 3.17: tt raw yields versus mass squared distribution. The histograms are our 
data. The curves are Gaussian fits. 
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Figure 3.18: raw yields versus mass squared distribution. The histograms are our 
data. The curves are Gaussian fits. 
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Figure 3.19: K raw yields versus mass squared distribution. The histograms are our 
data. The curves are Gaussian fits. 
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Figure 3.20: p raw yields versus mass squared distribution. The histograms are our 
data. The curves are Gaussian fits. 
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Figure 3.21: p raw yields versus mass squared distribution. The histograms are our 
data. The curves are Gaussian fits. 
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PT (GeV/c) 


d-l-Au Trigger 


0%-20% 


20%-40% 


40%- 100% 


P+P 


0.3-0.4 


2.861e -h 04 ± 169.4 


8922 ± 94.66 


8507 ± 92.38 


8519 ± 92.4 


1.715e + 04±131 


0.4-0.5 


2.139e -1- 04 ± 146.3 


6805 ± 82.49 


6458 ± 80.35 


6306 ± 79.41 


1.28e + 04 ± 113.1 


0.5-0.6 


1.611e -1-04 ±126.9 


5327 ± 72.98 


4873 ± 69.8 


4605 ± 67.86 


9189 ± 95.86 


0.6-0.7 


1.166e -1-04 ±108 


3831 ±61.9 


3550 ± 59.5 


3355 ± 57.92 


6362 ± 79.69 


0.7-0.8 


8447 ±91.91 


2837 ± 53.64 


2540 ± 50.4 


2387 ± 48.86 


4154 ± 64.5 


0.8-0.9 


5950 ±77.17 


2076 ±45.61 


1780 ± 42.2 


1646 ± 40.68 


2899 ± 53.87 


0.9-1 


4284 ± 65.46 


1446 ± 38.03 


1317 ± 36.31 


1171 ± 34.29 


1924 ± 44.01 


1-1.1 


3296 ± 57.47 


1123 ± 33.55 


1014 ±31.9 


897.5 ± 30.02 


1372 ± 37.77 


1.1-1.2 


2464 ± 49.88 


812.4 ±28.58 


762.4 ± 27.69 


650.8 ± 25.72 


1005 ±33.16 


1.2-1.4 


3136 ± 57.28 


1027 ± 32.54 


972.9 ±31.71 


828.8 ± 29.65 


1243 ± 39.78 


1.4-1.6 


1716 ± 45.79 


612.5 ±25.87 


539.3 ± 25.67 


422.7 ±24.22 


603.8 ± 30.49 


1.6-1.8 


1033 ± 39.74 


375.3 ± 21.21 


306.1 ± 19.59 


239.9 ± 48.55 


337.1 ± 28.82 



Table 3.4: tt raw signal table in minimum-bias, centrality selected d+Au collisions 
and minimum-bias p-|-p collisions. 



PT (GeV/c) 


d+Au Trigger 


0%-20% 


20%-40% 


40%-100% 


P+P 


0.4-0.5 


1410 ± 37.54 


417.2 ± 20.44 


420.9 ±20.52 


354.9 ± 18.84 


753.2 ± 27.44 


0.5-0.6 


1588 ± 39.85 


486.3 ± 22.06 


461 ±21.48 


435 ± 20.87 


729.3 ± 27 


0.6-0.7 


1499 ± 38.71 


465.9 ±21.59 


445.5 ±21.17 


395 ± 19.87 


710.1 ±26.65 


0.7-0.8 


1346 ± 36.69 


423.9 ± 20.59 


419.8 ±20.62 


335.7 ± 18.32 


579 ± 24.06 


0.8-0.9 


1105 ± 33.59 


369.7 ± 19.3 


317.9 ±18.43 


282.5 ± 17.03 


496.2 ± 22.44 


0.9-1 


969.1 ± 31.19 


283.9 ± 16.86 


305.1 ± 17.52 


258.7 ± 16.23 


381.2 ± 19.87 


1-1.1 


799.3 ± 28.41 


278.6 ± 16.79 


224 ± 15.04 


192.3 ± 14.04 


301.5 ± 18.39 


1.1-1.2 


656.7 ±26.24 


199.1 ± 14.33 


186.2 ± 14 


155.7 ± 12.94 


267.8 ± 18.71 


1.2-1.4 


1013 ± 34.43 


335 ± 19.51 


283.1 ± 17.71 


234.1 ± 17.68 


421.1 ±29.14 


1.4-1.6 


605.9 ±30.14 


191.1 ± 17.77 


174.1 ± 14.77 


148.8 ± 15.48 


238.8 ± 13.97 


1.6-1.8 


382.9 ± 28.9 











Table 3.5: K'^ raw signal table in minimum-bias, centrality selected d-l-Au collisions 
and minimum-bias p-|-p collisions. 



PT (GcV/c) 


d+Au Trigger 


0%-20% 


20%-40% 


40%-100% 


P+P 


0.4-0.5 


1341 ± 36.62 


411.6 ±20.29 


367.6 ± 19.19 


378 ± 19.44 


682.6 ±26.13 


0.5-0.6 


1498 ±38.7 


460 ±21.45 


411.5 ±20.32 


420.5 ± 20.51 


740.7 + 27.21 


0.6-0.7 


1410 ± 37.55 


436.2 ± 20.89 


398.4 ±19.97 


361.9 ± 19.02 


616.8 ±24.83 


0.7-0.8 


1207 ± 34.74 


366 ± 19.14 


349.7 ± 18.7 


350 ± 18.75 


557.4 ±23.61 


0.8-0.9 


1057 ± 32.66 


317.4 ± 18.12 


332.4 ± 18.37 


268.2 ± 16.66 


432.9 ± 20.93 


0.9-1 


863.7 ±29.42 


256.9 ± 16.09 


267.2 ± 16.43 


223.8 ± 15.03 


.368.9 ± 19.59 


1-1.1 


v,:i',.-2 ± 2r).:ij 


198.2 11.35 


183.1 13.01 


187.9 ± 13.88 


320.4 3z 19. 12 


1.1-1.2 


543 ± 23.92 


166.4 ± 13.14 


143 ± 12.22 


154.1 ± 12.89 


248.5 ± 18.84 


1.2-1.4 


895 ± 32.45 


302.1 ± 18.3 


258 ± 17.06 


206.7 ±16.17 


377.1 ± 26.67 


1.4-1.6 


645.5 ± 31.61 


202.4 ± 16.78 


141 ± 15.58 


166.6 ± 17.87 


237.6 ± 20.52 


1.6-1.8 


351.3 ±29.79 











Table 3.6: K raw signal table in minimum-bias, centrality selected d-l-Au collisions 
and minimum-bias p-|-p collisions. 
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PT (GeV/c) 


d-|-Au Trigger 


0%-20% 


20%-40% 


40%- 100% 


P+P 


0.4-0.5 


1377 ± 107.5 


412.6 ±40.73 


403.5 ±40.11 


422.7 ±41.42 


657.5 ±64.94 


0.5-0.6 


1527 ±89.71 


492.2 ±36.58 


428.6 ± 33.09 


424.8 ±32.89 


752.6 ±59.91 


0.6-0.7 


1156 ± 80.28 


437.9 ± 31.47 


120 ± 30.56 


101.3 ± 29.6 


701.9 ± 51.1 


0.7-0.8 


v.m> ± 73. (>8 


il0.() = 29. 13 


387.9 = 28.28 


385 ± 28.11 


(i7() ± r,5.3(> 


0.8-0.9 


1278 ± 72.45 


387.6 ± 28.57 


371.4 ±27.72 


312 ± 24.57 


498.8 ±45.17 


0.9-1 


1124 ± 68.87 


349.1 ±27.39 


322.1 ±25.87 


288.1 ±23.94 


448.3 ± 44.89 


1-1.1 


954.5 ± 62.39 


288.2 ± 24.38 


285.3 ± 24.23 


265.9 ± 23.11 


365 ±40.02 


1.1-1.2 


832.1 ± 58.34 


257.2 ± 23.06 


245.4 ±22.25 


213.3 ± 20.25 


273.8 ±34.08 


1.2-1.4 


1268 ± 72.57 


441.9 ±31.19 


362.6 ±27.1 


320.2 ± 24.77 


393.8 ±41.64 


1.4-1.6 


806.8 ±57.58 


306.9 ±26.38 


221.5 ±20.77 


190.8 ± 18.71 


210.3 ±31.59 


1.6-1.8 


540.8 ±23.27 


170.7 ± 13.06 


146.2 ± 12.14 


116.3 ± 10.81 


126 ± 11.31 


1.8-2 


314.2 ± 17.8 


119.2 ± 10.97 


81.7 ±9.764 


68.35 ± 9.093 


93.98 ± 10.15 


2-2.5 


388.1 ± 21.21 


148.4 ± 12.48 


135.7 ±12.01 


89.74 ± 10.02 


109 ± 12.33 


2.5-.'^ 


109.1 ± 12.92 


36.33 ± 6.809 


31.3 ± 8.488 


30.61 ± 7.187 


21.22 ± 5.122 


:-!- 1 


82.18 ± 12.30 











Table 3.7: p raw signal table in minimum-bias, centrality selected d+Au collisions 
and minimum-bias p-|-p collisions. 



PT (GcV/c) 


d+Au Trigger 


0%-20% 


20%-40% 


40%- 100% 


P+P 


0.4-0.5 


692.6 ±26.33 


215.1 ± 14.67 


183 ± 13.53 


202.8 ± 14.26 


421.7 ±20.56 


0.5-0.6 


1009 ± 31.76 


310.8 ± 17.63 


304 ± 17.43 


268.5 ± 16.39 


526.6 ±22.95 


0.6-0.7 


1098 ± 33.17 


317.9 ± 17.84 


305.6 ± 17.51 


327 ± 18.11 


561.5 ±23.71 


0.7-0.8 


1062 ± 32.-59 


340.2 ± 18.11 


307 ± 17.53 


285.5 ± 16.9 


135.1 ± 20.86 


0.8-0.9 


992.2 ± 31..") 


315.1 ^ 17.81 


281.1 ^ l(i.9() 


211. 1± 15.63 


370.1 ^ 19.1 


0.9-1 


827.5 ±28.76 


288.9 ± 17 


225 ± 15.01 


202.6 ± 14.24 


310.4 ± 17.62 


1-1.1 


724 ± 26.91 


240.2 ±15.5 


181.5 ±13.48 


192.2 ± 13.87 


246.4 ± 15.7 


1.1-1.2 


608.5 ± 24.67 


161.3 ±12.7 


184.3 ±13.61 


149.8 ± 12.25 


192 ± 13.87 


1.2-1.4 


914.9 ±30.24 


301.5 ± 17.36 


269.7 ±16.42 


214.1 ± 14.63 


269.6 ± 16.42 


1.4-1.6 


575.8 ± 24 


204.9 ± 14.32 


160.9 ±12.71 


120.5 ± 10.98 


138.6 ± 12.01 


1.6-1.8 


407.2 ± 20.18 


127.3 ± 11.29 


108.1 ± 10.43 


89.9 ±9.497 


100.6 ± 10.63 


1.8-2 


257.3 ± 16.26 


73.85 ±8.992 


92.22 ± 9.69 


46.81 ± 7.802 


71.23 ±8.92 


2-2.5 


305.6 ± 18.45 


114 ± 11.01 


83.43 ± 9.464 


77.84 ±9.16 


64.02 ± 10.44 


2.5-3 


111 ±12.79 


28.91 ± 6.26 


29.29 ± 8.869 


20.55 ±7.198 


25.71 ± 6.856 


3-4 


67.05 ± 11.87 











Table 3.8: p raw signal table in minimum-bias, centrality selected d-l-Au collisions 
and minimum-bias p-|-p collisions. 



Chapter 4 



Results 



4.1 TT^K^p and p spectra in d+Au and p+p colli- 
sions at mid-rapidity 
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Figure 4.1: The invariant yields of pions (filled circles), kaons (open squares), protons 
(filled triangles) and their anti-particles as a function of pr from d+Au and NSD p+p 
events at 200 GeV. The rapidity range was —0.5 < y < 0.0 with the direction of the 
outgoing Au ions as negative rapidity. Errors are statistical. 
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The invariant yields 2^~: fyj^^ of K^, p and p from both NSD p+p and 
minimum-bias d+Au events at mid-rapidity —0.5 < y < are shown in Figur jn| 
where is the corrected signal number per minimum-bias event in each p^ bin. 
N = '^''"Ivl^rxTal^rTtTl'tf.T"'' ^ where A^,,^ is the raw signal number in each pr bin, 
Ntotai is the total TOFr triggered events, factorl is the enhancement factor of TOFr 
trigger, factor2 is the TPC efficiency times TOFr matching efficiency, factorS is the 
background correction factor, factorA is the (Nch) bias factor, and factorb is the 
vertex efficiency times trigger efficiency and normalization factor. 

4.1.1 Systematic uncertainty 

For the invariant yield of vr^, K^, p and p, the average bin-to-bin systematic un- 
certainty was estimated to be of the order of 8%. The systematic uncertainty is 
dominated by the uncertainty in the detector response in Monte Carlo simulations 
(±7%). Additional factors contributing to the total systematic uncertainty include 
the background correction (±3%), the small rj acceptance of the TOFr (±2%), TOFr 
response (±2%), the correction for energy loss in the detector (^10 ± 10% at px < 0.6 
GeV/c for the p and p, much smaller for other species and negligible at higher p^), 
absorption of p in the material (±3%), and the momentum resolution correction 
(~ 5 ± 2%). The normalization uncertainties in d-l-Au minimum-bias and p+p NSD 
collisions are 10% and 14%, respectively [HI EE]- The charged pion yields are con- 
sistent with 71^ yields measured by the PHENIX collaboration in the overlapping px 
range fTlEZI. The invariant yields of vr^, K^, p and p in minimum-bias, centrality 
selected d+Au and minimum-bias p+p collisions, are listed in the tables in Appendix 
A with statistical errors and systematic uncertainties. 

4.2 Cronin effect 

Nuclear effects on hadron production in d+Au collisions are measured through com- 
parison to the p+p spectrum, scaled by the number of underlying nucleon-nucleon 
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Figure 4.2: The identified particle RdAu for minimum-bias and top 20% d+Au col- 
lisions. The filled triangles are for p + p, the filled circles are for 7r+ + tt^ and the 
open squares are for + K~ . Dashed lines are RdAu of inclusive charged hadrons 
from pS]. The open triangles and open circles are Rcp of p + p and 7r° in Au-|-Au 
collisions measured by PHENIX ^H] • Errors are statistical. The gray band represents 
the normalization uncertainty of 16%. 



inelastic collisions using the ratio 



RdAu 



£N/{2-npTdpTdy) 



TdAud^crfnei/ {'^T^PrdpTdy) ' 

where TdAu = {^bin) / o'inei describes the nuclear geometry, and d'^a^^^J {T.npTdpTdy) 
for p+p inelastic collisions is derived from the measured p+p NSD cross section. The 
difference between NSD and inelastic differential cross sections at mid-rapidity, as 
estimated from PYTHIA [TH], is 5% at low pt and negligible at pt > 1-0 GeV/c. 
Figure. 14.21 shows RdAu of 7r+ + tt^, + K~ and p + p for minimum-bias and 
central d+Au collisions. The systematic uncertainties on RdAu are of the order of 
16%, dominated by the uncertainty in normalization. The RdAu of the same particle 
species are similar between minimum-bias and top 20% d+Au collisions. In both 
cases, the RdAu of protons rise faster than RdAu of pions and kaons. We observe 
that the spectra of vr^, K"^, p and p are considerably harder in d-t-Au than those in 
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p+p collisions. The RdAu of the identified particles has characteristics of the Cronin 
effect j2Hl ISO] in particle production with RdAu less than unity at low pt and above 
unity at p^^^LO GeV/c. 

4.3 p + p/h ratio in d+Au and p+p collisions at 
middle pseudo-rapidity 
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Figure 4.3: Minimum-bias ratios of (p+p) over charged hadrons at — 0.5<?7<0.0 from 
y/sj^ = 200 GeV p+p (open diamonds), d+Au (filled triangles) and ^/s^ = 130 
GeV Au+Au jlHI (asterisks) collisions. Results of p + p collisions at ^/s^ = 1.8 
TeV pQ] are shown as open stars. Dashed lines are results of p/h^ ratios from 
■y/s^ = 23.8 GeV p+p (short-dashed lines) and p+W (dot-dashed) collisions [28j. 
Errors are statistical. 

Figure I^Sl depicts {p + p)/h, the ratio of p + p over inclusive charged hadrons as a 
function of px in d+Au and p+p minimum-bias collisions at y/s^ = 200 GeV, and 
p/h^ ratios in p+p and p+W minimum-bias collisions at y/s^ = 23.8 GeV [28 . Al- 
though the relative yields of particles and anti-particles are very different at ^/s < 40 
GeV due to the valence quark effects from target and projectile, the Cronin effects 
are similar. The systematic uncertainties on these ratios were estimated to be of the 
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order of 10% for p^^LO GeV/c, decreasing to 3% at higher px- At RHIC energies, 
the anti-particle to particle ratios approach unity {p/p = 0.81 ± 0.02 ± 0.04 in d+Au 
minimum-bias collisions) and their nuclear modification factors are similar. The dif- 
ference between RdAu at = 200 GeV forp + p and h can be obtained from the 
{p + p)/h ratios in d+Au and p+p collisions. Table Wl\ shows RdAu/ ^dAu determined 
by averaging over the bins within 1.2 < px < 3.0 GeV/c. At lower energy, the a pa- 
rameter in the power law dependence on target atomic weight A" of identified particle 
production falls with y/s [28^ . From the ratios of RdAu between p+p and h, we may fur- 
ther derive the ap-a^ for 1.2 < pr < 3.0 GeV/c to be 0.041±0.010(stat)±0.006(syst) 
under the assumptions that ok — ol-k and that (p+p)/7r and Kj-n are between 0.1 and 
0.4 in p+p collisions. This result is significantly smaller than the value 0.095 ± 0.004 
in the same pt range found at lower energies j28j . 

Table 4.1: {Nbin) from a Glauber model calculation, {p + p)/h averaged over the bins 
within 1.2 < pr < 2.0 GeV/c (left column) and within 2.0 < pr < 3.0 GeV/c (right 
column) and the RdAu ratios between p + p and h averaged over 1.2 < pt < 3.0 
GeV/c for minimum-bias, centrality selected d+Au collisions and minimum-bias p+p 
collisions. A p+p inelastic cross section of ainei = 42 mb was used in the calculation. 
For RdAu ratios, only statistical errors are shown and the systematic uncertainties are 
0.03 for all centrality bins. 



centrality 


{Nbin) 


{p + p)/h 


pP+P / D/l 

^dAu/ -"-dAu 


min. bias 
0-20% 
20-40% 
40 — 100% 


7.5 + 0.4 
15.0 + 1.1 
10.2 + 1.0 

^•'-'-0.3 


0.21 + 0.01 
0.21 + 0.01 
0.20 + 0.01 
0.20 + 0.01 


0.24 + 0.01 
0.24 + 0.02 
0.24 + 0.02 
0.23 + 0.02 


1.19 + 0.05 
1.18 + 0.06 
1.16 + 0.06 
1.13 + 0.06 


p+p 


1.0 


0.17 + 0.01 


0.21 + 0.02 





Also shown is {p+p)/h ratio from the Au+Au minimum-bias collisions at ^/s^ = 
130 GeV [18j. The {p + p)/h ratio from minimum-bias Au+Au collisions jTHj at a 
similar energy is about a factor of 2 higher than that in d+Au and p+p collisions for 
Pt'^2.0 GeV/c. This enhancement is most likely due to final-state effects in Au+Au 
collisions jHI 1211 113 1201 1221 12S| • The ratios show little centrality dependence in d+Au 
collisions, as shown in Table WTH and Figure For p^ < 2.0 GeV/c, the ratio in 



63 



0.5 
[0.45 

0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




dAu Minblas (/pp 1.19 ±0.05) 
dAu 0-20% (/pp 1.18 ±0.06) 
dAu 20-40% (/pp 1.16 ±0.06) 
dAu 40-100% (/pp 1.13 ±0.06) 



0.5 



1.5 



2.5 3 3.5 4 
Pt (GeV/c) 



Figure 4.4: Minimum-bias ratios of {p + p) over charged hadrons at —0.5 Kt] < 0.0 
from ^/s^ = 200 GeV minimum-bias and centrality selected d-l-Au collisions. Errors 
are statistical. 



p + p collisions at y/s^ = 1.8 TeV [8Q] is very similar to those in d+Au and p+p 
collisions at y/s^ = 200 GeV. 



4.4 K/n, p/n and anti-particle to particle ratios 



Centrality Bin 




X'^/ndf 


K-/K+ 


X'-'/ndf 


p/p 


X'^/ndf 


d+Au M.B. 


1.01 ±0.01 


0.88 


0.94 ± 0.02 


1.78 


0.81 ±0.02 


0.85 


0%-20% 


1.01 ±0.01 


0.80 


0.93 ± 0.03 


1.43 


0.80 ± 0.03 


0.70 


20%-40% 


1.00 ±0.01 


0.98 


0.91 ±0.03 


1.19 


0.79 ± 0.03 


1.14 


40%-100% 


1.02 ±0.01 


0.81 


1.02 ±0.03 


0.45 


0.78 ± 0.03 


0.70 


P+P 


1.00 ±0.01 


1.24 


0.98 ± 0.02 


0.71 


0.79 ± 0.03 


0.73 



Table 4.2: /n^, / and p/p ratios in p+p and d+Au minimum-bias collisions. 
Also shows in the table are the ratios in centrality selected d+Au collisions. Errors 
are statistical. 

Figure 031 shows the 7r^/7r+, / and p/p ratios as a function of Pt in d+Au 
and p+p minimum-bias collisions. It shows the anti-particle to particle ratios are fiat 
with Pt- The zero order polynominal function was used to fit the data and get the 
anti-particle to particle ratios. The results are list in Table 14.21 In centrality selected 
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Figure 4.5: 7r^/7r+, / and p/p ratios as a function of Pt in d+Au and p+p 
minimum-bias collisions. The open symbols are for p+p collisions and the solid 
symbols for d+Au collisions. Errors are statistical. 

d+Au collisions, the anti-particle to particle ratios are also fiat with px and show 
little centrality dependence. The results are also shown in the Table 14.21 
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Figure 4.6: i^/vr ratios as a function of pr in d+Au and p+p minimum-bias collisions. 
The open symbols are for p+p collisions and the solid symbols for d+Au collisions. 
Errors are statistical. 

The K/n andp/TT ratios are shown in Figure and Figure Wl\ individually. From 
the plots, the K/tt ratios increase with p^ in both d+Au and p+p collisions and the 
increasing trend is the same within our errors. The p/vr ratios increase with px in 
both d+Au and p+p collisions and the increasing in d+Au collisions is faster than 
that in p+p collisions. The trends of the K/n and p/vr as a function of pt show little 
centrality dependence in d+Au collisions. 
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Figure 4.7: p{P)/t^ ratios as a function of pt in d+Au and p+p minimum-bias col- 
lisions. The open symbols are for p+p collisions and the solid symbols for d+Au 
collisions. Errors are statistical. 

4.5 dN/dy, {Pt)^ cind model fits 

The spectra in minimum-bias and centrality selected d+Au collisions and also in p+p 
collisions are shown in Figure Figure and Figure 1^.1 Ul The spectra show little 
centrality dependence for each particle in d+Au collisions but harder than those in 
p+p collisions. The power law function was used to fit the spectra and get the dN/ dy 
and {pt)- The power law fit function is: 



1 d^N 



a(l + 



Pt 



2'npTdydpT ' ' (pr)^ 



(4.1) 



Where the parameter a is a constant value proportional to the mid-rapidity yield 
dN/dy, the parameter n is the order of the power law and (pt) is the mean value of 
the transverse momentum which is extracted from the fit. Figure E.lll shows power 
law fit to the spectra of minimum-bias d+Au and p+p collisions. Figure 14.121 shows 
power law fit to the spectra of 3 centrality selected d+Au collisions. The power law 
fit results are listed in Table 14.41 and Table 14.31 individually. The thermal model [HI] 
was also used to fit the spectra. The final dN/dy and (pt) are shown in Table and 
Table 14.51 respectively, which were obtained by averaging the results from the power 
law fit and thermal fit. Half of the differences in them are taken as the systematic 
errors due to the extrapolation to low p^ region. The errors in this table include the 
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Figure 4.8: The re-scaled vr"*" and vr spectra in minimum-bias, centrality selected 
d+Au collisions and also in p+p collisions. The errors are statistical. 



systematic uncertainties and statistical errors. 



Centrality Bin 


TV 


n+ 


K- 


K+ 


P 


P 


d+Au M.B. 


0.403 ± 0.004 


0.405 ± 0.004 


0.609 ± 0.009 


0.629 ± 0.009 


0.714 ±0.008 


0.677 ±0.010 


0%-20% 


0.421 ± 0.004 


0.421 ± 0.004 


0.626 ±0.018 


0.658 ± 0.016 


0.727 ±0.013 


0.705 ±0.014 


20%-40% 


0.408 ± 0.004 


0.411 ±0.004 


0.604 ±0.015 


0.625 ± 0.015 


0.725 ±0.013 


0.691 ±0.015 


40%-100% 


0.387 ±0.005 


0.391 ± 0.004 


0.589 ±0.016 


0.616 ± 0.016 


0.667 ±0.013 


0.646 ±0.014 


P+P 


0.357 ±0.004 


0.361 ±0.004 


0.571 ±0.013 


0.571 ± 0.013 


0.567 ±0.010 


0.569 ±0.012 



Table 4.3: (px) of tt^, 7r+, K~ , K^, p and p from power law fit in minimum-bias, 
centrality selected d+Au collisions and also in p+p collisions. The errors are from 
the power law fit. The unit of px is GeV/c. 



4.6 System comparison 

Figure E.13I shows the (px) of vr", K~ and p as a function of charged particle multi- 
plicity at mid-rapidity. From p+p to d+Au collisions, the (px) increase with charged 
particle multiplicity smoothly. We observed the {px) in 0%-20% d+Au collisions are 
larger than those in peripheral Au+Au collisions. The /tc^ and p/tc^ as a function 
of charged particle multiplicity at mid-rapidity are shown in Figure 1^.141 The /n^ 
and p/tt~ ratios were derived by taking the ratios of the dN/dy of K~ or p over the 
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Figure 4.9: The re-scaled and K spectra in minimum-bias, centrality selected 
d+Au collisions and also in p+p collisions. The errors are statistical. 



Centrality Bin 


7r 


77 + 


K- 


K+ 


P 


P 


d+Au M.B. 


5.078 ± 0.080 


5.032 ± 0.080 


0.685 ±0.013 


0.703 ±0.012 


0.466 ± 0.009 


0.594 ±0.019 


0%-20% 


10.657 ±0.190 


10.521 ±0.187 


1.448 ±0.085 


1.453 ±0.035 


0.972 ±0.026 


1.222 ±0.045 


20%-40% 


7.631 ±0.148 


7.515 ± 0.139 


0.988 ± 0.028 


1.051 ±0.027 


0.651 ±0.018 


0.842 ± 0.033 


40%-100% 


3.153 ±0.069 


3.024 ±0.060 


0.399 ±0.012 


0.379 ±0.011 


0.261 ±0.008 


0.338 ±0.014 


P+P 


1.524 ±0.027 


1.504 ±0.027 


0.166 ±0.004 


0.173 ±0.009 


0.113 ±0.003 


0.137 ±0.007 



Table 4.4: dN/dy of vr^, 7r+, ^ K^, p and p from power law fit in minimum-bias, 
centrality selected d+Au collisions and also in p+p collisions. The errors are from 
the power law fit. 



dN/dy of 7T~ in table ESI These ratios increase with charged particle multiplicity 
from p+p, d+Au to Au+Au collisions smoothly. The kinetic freeze out temperature 
Tkin and flow velocity {f3) from thermal fit as a function of charged particle multiplic- 
ity are shown in Figure 14.151 We can see the Tkin is flat from p+p to d+Au and then 
decreases from d+Au to Au+Au collisions and the increases from p+p, d+Au to 
Au+Au collisions. 
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Ccntrality Bin 




7r+ 


K- 


K+ 


P 


P 


d+Au M.B. 


0.420 ±0.019 


0.422 ± 0.019 


0.613 ± 0.025 


0.625 ± 0.025 


0.761 ± 0.056 


0.739 ± 0.069 


0%-20% 


0.435 ± 0.017 


0.436 ± 0.017 


0.627 ± 0.025 


0.646 ± 0.028 


0.774 ± 0.056 


0.761 ± 0.063 


20%-40% 


0.425 ± 0.019 


0.427 ±0.018 


0.610 ± 0.025 


0.622 ± 0.025 


0.766 ± 0.052 


0.744 ± 0.061 


40%- 100% 


0.405 ± 0.020 


0.408 ± 0.019 


0.591 ± 0.024 


0.608 ± 0.026 


0.715 ± 0.056 


0.703 ± 0.063 


P+P 


0.377 ±0.021 


0.379 ± 0.020 


0.565 ±0.023 


0.565 ± 0.023 


0.627 ±0.065 


0.634 ± 0.070 



Table 4.5: The final (pr) of tt , tt^, K , K^, p and p in minimum-bias, centrality 
selected d+Au collisions and also in p+p collisions. The unit of is GeV/c. 



Ccntrality Bin 


TV 


7r+ 


K- 


K+ 


P 


P 


d+Au M.B. 


4.731 ± 0.359 


4.668 ±0.356 


0.662 ± 0.040 


0.684 ± 0.039 


0.425 ± 0.044 


0.531 ± 0.067 


0%-20% 


10.063 ±0.628 


9.932 ±0.621 


1.383 ±0.095 


1.418 ±0.079 


0.896 ±0.084 


1.114 ±0.117 


20%-40% 


7.137 ±0.514 


7.074 ± 0.464 


0.952 ±0.059 


1.020 ±0.059 


0.603 ± 0.054 


0.765 ± 0.083 


40%- 100% 


2.925 ± 0.236 


2.829 ± 0.203 


0.387 ±0.023 


0.371 ±0.020 


0.238 ±0.025 


0.304 ±0.036 


P+P 


1.411 ±0.116 


1.400 ±0.108 


0.163 ±0.009 


0.168 ±0.010 


0.099 ±0.015 


0.120 ±0.018 



Table 4.6: The final dN/dy of tt , tt"*", K , K'^, p and p in minimum-bias, centrality 
selected d+Au collisions and also in p+p collisions. 
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Figure 4.10: The re-scaled p and p spectra in minimum-bias, centrality selected d+Au 
collisions and also in p+p coUisions. The errors are statistical. 
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Figure 4.11: The spectra of tt , tt"*", , , p and p in d+Au and p+p minimum-bias 
collisions. The curves are from power law fit. 
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Figure 4.12: The spectra of tt , tt"*", K , K^, p and p in three centrahty selected 
d+Au collisions. The curves are from power law fit. 
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Figure 4.13: (pt) as a function of dN/drj. The squared, circled and triangled symbols 
are from in p+p and Au+Au collisions. The cross, star and diamond are our 
data points in p+p and d+Au collisions. Statistic errors and systematic uncertainties 
have been added in quadrature. 
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Figure 4.14: K~ /t:~ and p/vr" as a function of dN/drj. The circled and triangled 
symbols are from [33] in p+p and Au+Au collisions. The star and diamond are our 
data points in p+p and d+Au collisions. Statistic errors and systematic uncertainties 
have been added in quadrature. 
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Figure 4.15: The kinetic freeze out temperature Tkin (left) and flow velocity (right) 
from thermal fit as a function of charged particle multiplicity. The circled symbols 
are from in p+p and Au+Au collisions. The star are our data points in p+p 
and d+Au collisions. Errors are systematic. 



Chapter 5 



Discussion 



5.1 Cronin effect 

The identified particle spectra in d+Au and p+p collisions not only provide the ref- 
erence for those in Au+Au collisions at 200 GeV, but also provide a chance to see 
the mechanism of the Cronin effect itself clearly. Cronin effect was observed 30 years 
ago [2H1- It is the enhancement of particle production at high pt- The enhancement 
was explained by initial multiple parton scattering. Also the recent experimental 
results of Cronin effect on inclusive charged hadron are consistent with the predic- 
tions based on initial multiple parton scattering jSO]. It suggests the suppression at 
intermediate px in Au+Au collisions is due to final state effects. However, the initial 
multiple parton scattering with the independent fragmentation function will result in 
the same Cronin effect for p{p) and for pions, while experimentally the Cronin effect 
for p{p) is larger than that for vr. That's to say the initial multiple scattering with 
the independent fragmentation function can't account for the Cronin effect observed. 
Maybe in the initial multiple parton scattering, the broadening for gluon and for 
quark/antiquark are not the same |H8j. Or maybe the fragmentation processes in 
p+A collisions are not the same as those in p+p collisions [21]. Whether the Cronin 
effect is initial state effect or final state effect will be discussed below. 
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5.1.1 Model comparison: initial state effect? 

The initial multiple parton scattering model predicts that the Cronin effect on deuteron 
beam outgoing side is larger than that on Au beam outgoing side since the deuteron 
traverses a much larger nucleus ^7\. Figure I^TTl (left) shows the predictions for the 
Cronin effect at different rapidity range. The different curves correspond to the pre- 
diction results from different shadowings. The ?/ = 1 is on the deuteron beam outgoing 
side. The ?/ = — 1 is on the Au beam outgoing side. The ?/ = is at mid-rapidity. We 
can see that the RdAu on deuteron beam side {y = 1) increases faster than that on 
Au beam side {y = —1). If we take the ratio of RdAu on Au beam side over RdAu on 
deuteron beam side, it will result in a minimum value at ~ 3.5 GeV/c, as shown 
in the curves on the right plot of Figure 15.11 The solid symbol on the right plot of 
Figure EH) represents the data points jHO], which is the ratio of RdAu on Au beam side 
at — 1 < ?7 < —0.5 over RdAu on deuteron beam side at 0.5 <!]<!. We observe the 
1] asymmetry from experiment reaches a maximum value firstly and then decreases. 
This is different from the predictions. That means, the model based on initial mul- 
tiple parton scattering only, can't reproduce the experimental results. Recently, Qiu 
and Vitev have come up with the idea of coherent multiple scattering and applied it 
to the RHIC experiments |ST] . In this picture, the hard probe may interact coherently 
with many low x parton inside different nucleons inside the nucleus. As a result, this 
process will lead to the suppression of the total cross section. This coherent effect 
wiU play an important role in p+A collisions at forward rapidity. In the deuteron 
outgoing beam direction, the coherent effect is non-negligible since the Au nucleus 
is big while on the Au side, the coherent effect is not big since the deuteron is of a 
small size. This will result in bigger suppression on the deuteron side than on the Au 
side. It may qualitatively reproduce the data. This coherent multiple scattering is a 
final-state effect. 

As we all known, in Au+Au collisions, the suppression at intermediate pt can be 
reproduced by the initial multiple scattering and jet quenching qualitatively [Hlini- 
However, the model based on the initial multiple scattering, jet quenching and inde- 
pendent fragmentation will result in the same suppression for baryons and mesons at 
intermediate px in Au+Au collisions. Experimentally Rep for baryons are larger than 
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Figure 5.1: (left) The Cronin effect at different rapidity as a function of pt- Tlie 
different curve in each panel shows the different shadowing. This figure is from |87]. 
(right) The 7] asymmetry of the Cronin effect: the ratio of Cronin effect in Au beam 
outgoing direction over the Cronin effect in deuteron beam outgoing direction. This 
figure is from pO] . 



Rep for mesons at intermediate pt- This difference can be reproduced by coalescence 



or recombination models j^U 1221 122] • Recently the recombination model [121 has been 
applied to d+Au system to see whether it can reproduce the Cronin effect or not. 

With the help of Prof. C.B. Yang [93], I also compare our pion and proton 
spectra in d+Au collisions with the recombination model [22! • ■'■^ ^^e following the 
recombination model [22] will be discussed and the comparison between the data and 
the model will be presented in detail. 



5.1.2 Model comparison: recombination 

The inclusive distribution for the production of pions can be written in the recombi- 
nation model jn2l! when mass effects are negligible, in the invariant form 



p- 



dp 



dpi dp2 
Pi P2 



FggiPl,P2)RAPl,P2,p), 



(5.1) 
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where Fgg{pi,p2) is the joint distribution of a g and q at pi and p2, and Rt,{j>i,P2,p) is 
the recombination function for forming a pion at p: Rt,{j)i,P2, p) = {piP2/p)S{Pi +P2 — 
p). Fqq depends on the coUiding hadron/nuclei. In general, Fqq has four contributing 
components represented schematically by 

Fqq = TT + TS + {SS)i + {SS)2 (5.2) 

where T denotes thermal distribution and S shower distribution. {SS)i signifies 
two shower partons in the same hard-parton jet, while {88)2 stands for two shower 
partons from two nearby jets |92j . 

Yoi p + A collisions it may not be appropriate to refer to any partons as thermal in 
the sense of a hot plasma as in heavy-ion collisions. Here in d+Au collisions, the sym- 
bol T represents the soft parton distribution at low kx- At low px the observed pion 
distribution is exponential; we identify it with the contribution of the TT term |i92j . 

where T is the inverse slope. We shall determine C and T by fitting the d+Au data 
at low pt- The pion spectra for different centralities can be calculated from thermal- 
thermal {pioutt), thermal-shower [piorits) and shower-shower {pioHss) contributions by 
using parameters C and Nbin. dN/prdpT = C x C x pionu + 2.5 x C x Nun x piorits + 
2.5 X Nun X piouss, where C is determined by fitting the d+Au data at 0.4 < px < 1-0 
GeV/c, and Nhm is the number of binary collisions. The data points of pioriu, piorits 
and pioriss are from Prof. CB. Yang [HSj. The C values for minimum-bias, 0-20%, 
20-40% and 40-~100% d+Au coUisions are 8.85, 13.08, 10.96 and 6.84 individually. 
The T value of 0.21 GeV is used in the low px fit. Figure shows the tt^ spectra in 
d+Au collisions as well as those from recombination model. This figure shows that 
the recombination model can reproduce the spectra of pion in minimum-bias and 
centrality selected d+Au collisions. 

The invariant inclusive distribution for proton formation at midrapidity in the 
recombination model [221 

^QdN^^ /"^^i?(p^^P2,P3)^p(Pl,P2,P3,P) (5.4) 

dp J pi P2 
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Figure 5.2: (left) The invariant yield for 7r+ at 0%-20% d+Au collisions as a function 
of Pt- The open circles are our data points. The curves are the calculation results 
from recombination model. Sum represents the total contribution from recombina- 
tion model. Thermal-thermal represents the soft contribution. The thermal- shower 
represents the contribution from the interplay between soft and hard components. 
The shower-shower represents the hard contribution, (right) The invariant yields for 
tt"^ in minimum-bias and centrality selected d-l-Au collisions as a function of pt- The 
symbols represent our data points. The curves on the top of the symbols are the 
corresponding calculation results from recombination model. 

where all momentum variables pi and p are transverse momenta, and p^ denotes the 
energy of the proton. F{pi,p2,P3) is the joint distribution of u,u, and d quarks at 
Pi,P2 and p3, respectively. Rp{pi,P2,P3,p) is the recombination function for a proton 
with momentum p. We write schematically 

F = TTT + TTS + TSS + SSS (5.5) 

where all the shower partons S are from one hard parton jet. Shower partons from 
different jets are ignored here for RHIC energies. In d-l-Au coUisions, T denotes the 
soft partons that are not associated with the shower components of a hard parton. 
The SSS term is regarded as the fragmentation of a hard parton into a proton. The 
TTT term comes entirely from the soft partons, while TTS and TSS accounts for 
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Figure 5.3: (left) The invariant yield for p at 0%-20% d+Au collisions as a function 
of pt- The open circles are our data points. The curves are the calculation results 
from recombination model. Sum represents the total contribution from recombination 
model. TTT represents the soft contribution. The TTS+TSS represents the contri- 
bution from the interplay between soft and hard components. The SSS represents the 
hard contribution, (right) The invariant yields for p in minimum-bias and centrality 
selected d+Au collisions as a function of pt- The symbols represent our data points. 
The curves on the top of the symbols are the corresponding calculation results from 
recombination model. 

the interplay between the soft and shower partons. The soft contribution to the 
proton spectrum arising from TTT recombination is 

pdp 6 5(a + l,7 + l)5(a + l,a + 7 + 2) ^'^ 

. Where C and T are determined by fitting the proton spectra at low pr, the a is equal 
to 1.75, 7 is equal to 1.05, B{x,y) is the beta function [92j. For the invariant yield of 
proton, there are 4 different contributions: soft-soft-soft (protonut), soft-soft-shower 
{protoritts), soft-shower-shower [protoritss), and shower-shower-shower [protongss)- The 
total contributions are dN/pxdpT = C x C x C x protonut + C xC x N^n x protonus + 
C X Nbin X protoritss + Nbm X protousss, where C is determined by fitting the d+Au 
data at 0.5 < pt < 1-5 GeV/c, and Ni,in is the number of binary collisions. The data 
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points of protorittt, protoutts, protoutss and protoUsss are from Prof. C.B. Yang PHj. 
The C values for minimum-bias, 0-20%, 20-40% and 40-~100% d-l-Au collisions are 
9.67, 12.34, 10.92 and 7.91 individually. The T value of 0.21 GeV is used in the low 
Pt fit. Figure EISl shows the proton spectra in d-|-Au collisions as well as those from 
recombination model. This figure shows that the recombination model can reproduce 
the spectra of proton in minimum-bias and centrality selected d-|-Au collisions. From 
the comparison between our data and the calculation results from the recombination 
model, we know that the recombination model actually can reproduce both the proton 
and pion spectra in d-l-Au collisions, while as we have mentioned above, the initial 
multiple parton scattering model jSOl with independent fragmentation can't reproduce 
the difference of Cronin effect between proton and pion. Besides, the initial multiple 
parton scattering model with independent fragmentation can't reproduce the t] asym- 
metry of the Cronin effect. In the recombination model jH^! the number of such soft 
partons on the Au outgoing side is larger than that on the deuteron outgoing side. 
This will result in the Cronin effect on the Au side larger than that on the deuteron 
side jn2]- Qualitatively the recombination model can reproduce the r] asymmetry of 
the Cronin effect. As we know, the recombination model is a final-state effect model. 
These all seem to indicate that the Cronin effect is not initial-state effect only. The 
final-state effect plays an important role too. To directly confirm the Cronin effect 
is initial or final state effect, it's necessary for us to compare the Cronin effect of 
Drell-Yan process with those of pion, kaon and proton. I will come to this later. 

5.1.3 Integral yield RdAu- shadowing effect? 

The initial multiple elastic scattering only changes the pt distribution while the total 
cross section should not change. Thus we can look at the integral yield dN/dy RdAu, 
which are measured through comparison to the integral yield dN/ dy in p+p collisions, 
scaled by the number of binary collisions Nbin- Figure ED shows that integral yield 
RdAu of pion, kaon and proton as a function of dN/ drj in minimum-bias and centrality 
selected d-l-Au collisions at mid-rapidity. The integral yield RdAu of pion and kaon 
are less than 1 while that of proton is close to 1. This may be the indication of 
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Figure 5.4: Integral yield RdAu as a function of dN/dt] in minimum-bias and centrality 
selected d+Au collisions at mid-rapidity. Statistic errors and systematic uncertain- 
ties have been added in quadrature. The shadowing represents the normalization 
uncertainty. 

shadowing effect at 200 GeV. The integral yield RdAu for proton is larger than that 
for kaon and a little bit more larger than that for pion. This may be the indication 
that the shadowing effect is mass dependent at 200 GeV. 



5.1.4 Initial or final state effect: Drell-Yan process 

In order to see the Cronin effect is initial or final state effect, we may look into the 
Drell-Yan process since there is little final state effect in Drell-Yan process. If there 
is no enhancement at high pt for Drell-Yan process, the enhancement for iT,K,p is 
due to final-state effect. Figure E31 shows the integral yield Cronin ratio as a function 
of atomic weight at p-A fixed target experiment [2^. The proton incident energy is 
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800 GeV. We can see there is no enhancement for Drell-Yan process. However, this 
is the total cross section while what we want to compare is Cronin ratio as a function 
of pt- It will be better if we have the pt dependence of Cronin ratio of Drell-Yan 
process. However, at the same pt range with the same proton incident energy, the 
Cronin ratio of Drell-Yan is not available in p+A collisions. It's hard to compare the 
Cronin ratio of Drell-Yan process with those of vr, K,p. 
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Figure 5.5: Integral yield Cronin ratio as a function of atomic weight in p+A fixed 
target experiment for Drell-Yan process, etc. This figure is from 95]. 



5.2 Baryon excess in Au+Au collisions 

Now let's come to another important physics from d+Au collisions. We know that 
the [p + p)/h ratio from minimum-bias Au+Au collisions JH] at a similar energy is 
about a factor of 2 higher than that in d+Au and p+p collisions for pt^2.0 GeV/c. 
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This enhancement is most hkely due to final-state effects in Au+Au colhsions. There 
are many models trying to explain this baryon excess in Au+Au collisions jHl 122 ^1 
1201 1221 123 • In the following baryon production mechanism will be discussed. 

5.2.1 p/p ratio vs pt 




- Au+Au M.B. from PHENIX 
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Figure 5.6: p/p ratio as a function of pr in d+Au and p+p minimum-bias colli- 
sions. The open squared symbols are for p+p collisions and the solid circled sym- 
bols for d+Au collisions. The triangled symbols represent the result from Au+Au 
minimum-bias collision 14 . The curve is the pQCD calculation results from [21] in 
p+p collisions. Errors are statistical. 

In 200 GeV Au+Au collisions, p/p ratio was observed to be flat with px till 
intermediate pt range [I3|, as shown in Figure ESI The baryon junction model [21] 
tried to explain it by using junction anti-junction production with jet quenching, 
on the basis of pQCD calculation [21] where the p/p ratio decreases with pj- in p+p 
collisions. The curve from pQCD calculation [2^ is also shown in Figure|SISl However, 
p/p ratios in d+Au and p+p collisions in our data show to be flat with px within 
errors. Anyway, the precise measurement with more statistics in p+p and d+Au 
collisions is needed to address this issue. 
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5.2.2 Baryon production at RHIC: multi-gluon dynamics? 
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Figure 5.7: Minimum-bias ratios of {p + p) over charged hadrons at —0.5 Ki] < 0.0 
from y/s^ = 200 GeV p+p (open diamonds) and d+Au (filled triangles) collisions. 
Also shown are the {p + p)/h ratios in e+e" collisions at ARGUS PHI- The solid line 
represents the {p + p)/h ratio from three gluon hadronization while the dashed line 
for the ratio from quark and antiquark fragmentation |S6j. Errors are statistical. 

Let's compare the {p + p)/h ratio in p+p, d+Au and Au+Au collisions at RHIC 
energy 200 GeV with the ratio in e^e~ collisions at ARGUS ^H]- Using the ARGUS 
detector at the e^e~ storage ring DORIS II, the inclusive production of pion, kaon and 
proton in multihadron events at 9.98 GeV and in direct decays of the T(IS') meson 
were investigated jHni- Multihadron final states in e^e~ annihilation are produced 
via quark and antiquark fragmentation, and those from direct T(IS') decays originate 
from the hadronization of three gluons [HH] • Figure 15.71 shows the {p + p)/h ratio 
in 200 GeV p+p collisions together with the ratio in e~^e~ collisions at ARGUS [96j. 
The plot shows that the {p + p)/h ratio from three gluon hadronization is a factor of 3 
higher than that from quark and antiquark fragmentation at ARGUS. Our data from 
200 GeV p+p collisions is close to {p + p)/h ratio from 3 gluon hadronization. This 
may be the indication that in the heavy ion collisions at RHIC energy, multi-gluon 
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hadronization plays an important role for the particle production. 



Chapter 6 

Conclusion and Outlook 



6.1 Conclusion 

In summary, we have reported the identified particle spectra of pions, kaons, protons 
and anti-protons at mid-rapidity from 200 GeV minimum-bias, centrality selected 
d+Au collisions and NSD p+p collisions. The time-of-flight detector, based on novel 
multi-gap resistive plate chamber technology, was used for particle identification. This 
is the first time that MRPC detector was installed to take data as a time-of-fiight 
detector in the collider experiment. The calibration method was set up in the STAR 
experiment for the first time and has been applied to the data taken later successfully. 
The intrinsic timing resolution of the MRPC was 85 ps after the calibration. In 2003 
run, the pion/kaon can be separated up to transverse momentum 1.6 GeV/c while 
proton can be identified up to 3.0 GeV/c. 

The spectra of tt^, K"^, p and p in d-l-Au and p-|-p collisions provide an important 
reference for those in Au-|-Au collisions. The initial state in d-l-Au collisions is similar 
to that in Au-|-Au coUisions, and, it's believed that the quark-gluon plasma doesn't 
exist in d-l-Au collisions. These results from d-l-Au collisions are very important for 
us to judge whether the quark-gluon plasma exists in Au-|-Au collisions or not and to 
understand the property of the dense matter created in Au-|-Au coUisions. We observe 
that the spectra of tt^, K^, p and p are considerably harder in d-l-Au than those in 
p-|-p coUisions. In ^/s]^ — 200 GeV d-l-Au coUisions, the RdAu of protons rise faster 
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than RdAu of pions and kaons. The RdAu of proton is larger than 1 at intermediate 
Pt while the proton production follows binary scaling at the same pt range in 200 
GeV Au+Au collisions. These results further prove that the suppression observed 
in Au+Au collisions at intermediate and high px is due to final state interactions 
in a dense and dissipative medium produced during the collision and not due to 
the initial state wave function of the Au nucleus. Additionally, the particle-species 
dependence of the Cronin effect is found to be significantly smaller than that from 
lower energy p+A collisions. In y/s^ = 200 GeV d+Au collisions, the ratio of the 
nuclear modification factor RdAu between {p + p) and charged hadrons (h) in the 
Pt range 1.2 < pt < 3.0 GeV/c was measured to be 1.19 ± 0.05(stat)±0.03(syst) 
in minimum-bias collisions. Both the RdAu values and {p + p)/h ratios show little 
centrality dependence, in contrast to previous measurements in Au+Au collisions at 
y/SNN = 130 and 200 GeV. The ratios of protons over charged hadrons in d+Au 
and p+p collisions are found to be about a factor of 2 lower than that from Au+Au 
collisions, indicating that the relative baryon enhancement observed in heavy ion 
collisions at RHIC is due to the final state effects in Au+Au collisions. 

The identified particle spectra in d+Au and p+p collisions not only provide the 
reference for those in Au+Au collisions, but also provide a chance to see the mecha- 
nism of the Cronin effect itself clearly. Usually the Cronin effect has been explained 
to be the initial state effect only since 1970s PU|. However, we compare our pion 
and proton spectra in minimum-bias and centrality-selected d+Au collisions with the 
recombination model jHll- The recombination model can reproduce both the pion 
spectra and proton spectra. This recombination model is built on the hadroniza- 
tion process, which is a final-state effect, while the initial multiple parton scattering 
model can't reproduce the difference of the Cronin effect between pions and pro- 
tons. From these comparisons, we conclude that the Cronin effect in ^,/s^ = 200 
GeV d+Au collisions is not the initial state effect only, and that final state effect 
plays an important role. 

The integral yield dN/dy and {pt) in p+p and d+Au collisions were estimated 
from the power law fit and thermal model fit. The integral yield RdAu of vr^, K^, p 
and p are observed to be smaller than 1 while those of p and p are close to 1. The 
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/-n^ , / and pjp ratios as a function of pt are observed to be flat with pt 
within the errors in d+Au and p+p minimum-bias colhsions and show httle centrahty 
dependence in d+Au colhsions. The integral yield ratios of K~ j-n" and p/vr~ as a 
function of dN/dr] were also presented in p+p and d+Au collisions. 

6.2 Outlook 

For the outlook, I will discuss whether the Cronin effect is mass dependent or baryon/meson 
dependent at 200 GeV. What other physics topic have we done from MRPC-TOFr in 
d+Au and p+p collisions in 2003 run? If we have the full time-of-fiight (FuU-TOF) 
coverage, what can we do? Also I will discuss a little bit about the low energy 63 
GeV Au+Au run. 

6.2.1 Cronin effect at 200 GeV: Mass dependent or baryon/meson 
dependent? 

We know that recombination model can reproduce the spectra of pions and protons in 
d+Au collisions. Also the Rcp of identified particles in Au+Au collisions suggest that 
the degree of suppression depends on particle species(baryon/meson) at intermediate 
Pt- Does the Cronin effect in 200 GeV d+Au collisions depend on the particle species 
(baryon/meson) or depend on the particle mass? From our data, it shows the Cronin 
effect for proton is bigger than those for pion and kaon. And the Cronin effect of 
pion shows little difference from that of kaon at pt < 1.5 GeV/c. In order to see the 
Cronin effect is baryon/meson dependent or mass dependent, we can compare the 
Cronin effect of proton with those of K* and (f) since the mass of K* and (p are close 
to that of proton while K* and are mesons and proton is a baryon. The preliminary 
results show that the Cronin effect of K* and PZ] are similar to that of pion and 
different from that of proton. However, the final results from K* and are needed 
to confirm this issue. 
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6.2.2 Electron PID from MRPC-TOFr 

The production and spectra of hadrons with heavy flavor are sensitive to initial con- 
ditions and the later stage dynamical evolution in high energy nuclear collisions, 
and may be less affected by the non-perturbative complication in theoretical calcu- 
lations |lUOj . Charm production has been proposed as a sensitive measurement of 
parton distribution function in nucleon and the nuclear shadowing effect by system- 
atically studying p+p, and p+A collisions |lUlj . The relatively reduced energy loss of 
heavy quark traversing a quark-gluon plasma will help us distinguish the medium in 
which the jet loses its energy |lU2j . A possible enhancement of charmonia (J/^E') pro- 
duction can be present at RHIC energies |lU3j due to the coalescence of the copiously 
produced charm quarks |99|. 

X 
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D 0.5 1 1^ 2 2.5 3 

Figure 6.1: dE/dx in TPC versus p without(the upper panel) or with (the lower 
pannel) the TOFr velocity cut — 1| < 0.03. The insert shows dE/dx distribution 
for 1 < p < 1.5 GeV/c. 

The recent STAR results on the absolute open charm cross section measurements 
from direct charmed hadron reconstruction |^ in d+Au collisions and electrons 
from charm semileptonic decay in both p+p and d+Au collisions at 200 GeV were 
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presented jHHl- Based on the capability of hadron identification [HH] from the MRPC- 
TOFr tray in 2003, electrons could be identified at low momentum {p^ < 3 GeV/c) 
by the combination of velocity from TOFr and the particle ionization energy 
loss (dE/dx) from TPC j38|. Figure [HH] shows that the electrons are clearly identified 
as a separate band in the dE/dx versus momentum {p) with a selection on j3 at 
1| < 0.03 in d+Au collisions. At higher (2-4 GeV/c), negative electrons were 
also identified directly by TPC since hadrons have lower dE/dx due to the relativisitic 
rise of electron dE/dx. Based on the clear electron identification, the open-charm- 
decayed electron spectra was derived Combined with measurement from 
TPC, the total charm cross section was obtained ^\ in d+Au collisions. 

6.2.3 Full-TOF Physics 

Based on the hadron PID and electron PID of MRPC-TOFr in 2003, we can imagine 
how many physics we can do if we have full time-of- flight coverage based on MRPC 
technology. The proposal [7^ for large area time-of-fiight system for STAR has been 
proposed. Since the pion/kaon can be separated up to transverse momentum 1.6 
GeV/c and proton can be identified up to 3.0 GeV/c from time-of-fiight system. The 
resonance spectra measured from hadronic decay will be extended to much higher px- 
The direct open charm spectra from its hadronic decay channel will reach higher pt 
with much more precise measurement. Since the electron can be clearly identified up 
to transverse momentum 3~4 GeV/c by the combination of velocity from TOFr 
and the particle ionization energy loss (dE/dx) from TPC jSHj, the electron 
spectra from charm-semi-leptonic decay will be measured precisely. As we know that 
the measurement of the di-leptonic decays of vector mesons are very difficult since 
the branch ratios are too small and it's really hard to subtract the background. But 
with TOF upgrading together with the SVT and micro-vertex detector upgrading, 
the di-leptonic decays of vector mesons will be measured much more easily, which 
will bring the direct information of QGP since the electron is a lepton and the cross 
section of interaction between electrons and hadrons is little. Thus we can see directly 
the property of quark-gluon plasma such as the temperature and the chiral symmetry 
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restoration. This will be the most interesting and meaningful thing for the QGP 
search jM] . Besides, there are many other physics topics [71j such as identified particle 
correlation and fluctuation, particle composition of jet fragmentation, and anti-nuclei 
etc. 

6.2.4 63 GeV Au+Au collisions at RHIC 

The bulk properties such as elliptic flow V2 and particle production show smooth 
trend from AGS, SPS to RHIC energy. One energy point y/s^ = 63 GeV, which is 
between SPS and full RHIC energy, was selected since high quality charged-particle 
and 7r° inclusive spectra have been measured in p+p collisions at 63 GeV at Inter- 
secting Storage Rings (ISR) and will serve as the reference spectra for computing the 
nuclear modification factor for Au+Au collisions measured at the same energy. The 
V2 and particle production at 63 GeV will be studied at RHIC. Besides, the nuclear 
modification factor as a function of pt will also be studied in this collision system 
in which the hard scattering component has been significantly reduced. The results 
from 63 GeV Au+Au collisions will be helpful for us to understand the property of 
dense medium created in 200 GeV Au+Au collisions. 



Appendix A 

Tables of the Invariant 

Yields 



Pt 


Pt width 


M.B. 


0%-20% 


3.50e-01 


l.OOe-01 


3.14e -h 00 ± 3.65e - 02 ± 2.51e - 01 


6.70e -h 00 ± 9.69e - 02 ± 5.36e - 01 


4.50e-01 


l.OOe-01 


1.82e -t- 00 ± 2.45e - 02 ± 1.45e - 01 


3.88e -h 00 ± 6.51e - 02 ± 3.10e - 01 


5.50e-01 


l.OOe-01 


1.08e -h 00 ± 1.73e - 02 ± 8.67e - 02 


2.38e -h 00 ± 4.68e - 02 ± 1.90e - 01 


6.50e-01 


l.OOe-01 


6.18e - 01 ± 1.22e - 02 ± 4.94e - 02 


1.38e -h 00 ± 3.28e - 02 ± l.lOe - 01 


7.50e-01 


l.OOe-01 


4.13e - 01 ± 9.15e - 03 ± 3.31e - 02 


9.53e - 01 ± 2.56e - 02 ± 7.63e - 02 


8.50C-01 


l.OOc-01 


2.66e - 01 ± 6.98e - 03 ± 2.13e - 02 


e.lle - 01 ± 1.93e - 02 ± 4.88e - 02 


9.50C-01 


l.OOc-01 


1.70e - 01 ± 5.36e - 03 ± 1.36e - 02 


3.79e - 01 ± 1.44e - 02 ± 3.03e - 02 


1.05C+00 


l.OOc-01 


1.14e - 01 ± 4.07e - 03 ± 9.09e - 03 


2.68e - 01 it 1.15e - 02 it 2.14e - 02 


1.15C+00 


l.OOc-01 


7.70e - 02 ± 3.27e - 03 ± 6.16e - 03 


1.71e - 01 di 8.82e - 03 ± 1.37e - 02 


1.30e+00 


2.00e-01 


4.43e - 02 ± 1.66e - 03 ± 3.55e - 03 


1.06e - 01 it 4.71e - 03 ± 8.45e - 03 


1.50e+00 


2.00e-01 


2.10e - 02 it 1.07e - 03 ± 1.68e - 03 


4.65e - 02 it 2.96e - 03 it 3.72e - 03 


1.70e+00 


2.00e-01 


1.05e - 02 ± 7.19e - 04 ± 8.40e - 04 


2.29e - 02 it 1.92e - 03 it 1.83e - 03 



Table A.l: tt^ spectra in minimum-bias and 0-20% d-f-Au collisions. The unit of pt 
and Pt width is GeV/ c. 
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PT 


Pt width 


20%-40% 


40%-100% 


3.50e-01 


l.OOe-01 


4.73e -h 00 ± 6.94e - 02 ± 3.78e - 01 


1.93e -h 00 ± 2.85e - 02 ± 1.55e - 01 


4.50e-01 


l.OOe-01 


2.80e -h 00 ± 4.73e - 02 ± 2.24e - 01 


l.lOe -h 00 ± 1.89e - 02 ± 8.81e - 02 


5.50e-01 


l.OOe-01 


1.68e -h 00 ± 3.36e - 02 ± 1.34e - 01 


6.43e - 01 ± 1.31e - 02 ± 5.14e - 02 


6.50e-01 


l.OOe-01 


9.50e - 01 ± 2.30e - 02 ± 7.60e - 02 


3.68e - 01 ± 9.03e - 03 ± 2.94e - 02 


7.50C-01 


l.OOc-01 


6.41e - 01 ± 1.76e - 02 ± 5.13e - 02 


2.40e - 01 ± 6.80e - 03 ± 1.92e - 02 


8.50C-01 


l.OOc-01 


4.20e - 01 ± 1.36e - 02 ± 3.36e - 02 


1.52e - 01 ± 5.09e - 03 ± 1.22e - 02 


9.50e-01 


l.OOc-01 


2.58e - 01 ± l.OOe - 02 ± 2.07e - 02 


l.Ole - 01 ± 3.97e - 03 ± 8.05e - 03 


1.05C+00 


l.OOc-01 


1.80e - 01 ± 7.92e - 03 ± 1.44e - 02 


e.lOe - 02 ± 2.84e - 03 ± 4.88e - 03 


1.15C+00 


l.OOc-01 


1.21e - 01 ± 6.35e - 03 ± 9.72e - 03 


4.38e - 02 ± 2.37e - 03 ± 3.50e - 03 


1.30(-+00 


2.00(--01 


6.19r - 02 ± 3.02k - 03 ± 5.19c - 03 


2.10c - 02 ± 1.16c - 03 ± 1.92c - 03 


1.5()c+()() 


2.0()c-()l 


3.28c - 02 ± 2.().5c - 1)3 = 2.(>3c - 03 


1.03c - 02 = 7.12c - 01 = 8.27c - 04 


1.70e+00 


2.00e-01 


1.57e - 02 ± 1.31e - 03 ± 1.26e - 03 


4.92e - 03 ± 4.49e - 04 ± 3.94e - 04 



Table A.2: 7r+ spectra in 20-40% and 40-100% d+Au collisions. The unit of pt and 
Pt width is GeV/c. 



PT 


Pt width 


P+P 


3.50C-01 


l.OOc-01 


9.71e - 01 ± 1.21e - 02 ± 7.76e - 02 


4.50C-01 


l.OOc-01 


5.32e - 01 ± 7.84e - 03 ± 4.26e - 02 


5.50C-01 


l.OOc-01 


3.14e - 01 ± 5.46e - 03 ± 2.51e - 02 


6.50e-01 


l.OOe-01 


1.74e - 01 ± 3.72e - 03 ± 1.40e - 02 


7.50e-01 


l.OOe-Ol 


1.08e - 01 ± 2.65e - 03 ± 8.64e - 03 


8.50e-01 


l.OOe-01 


6.42e - 02 ± 1.89e - 03 ± 5.14e - 03 


9.50e-01 


l.OOe-01 


4.03e - 02 ± 1.42e - 03 ± 3.22e - 03 


1.05e+00 


l.OOc-01 


2.40e - 02 ± 9.94e - 04 ± 1.92e - 03 


1.15e+00 


l.OOc-01 


1.55e - 02 ± 7.71e - 04 ± 1.24e - 03 


1.30C+00 


2.00C-01 


8.19e - 03 ± 3.86e - 04 ± 6.55e - 04 


1.50e-|-00 


2.00e-01 


3.77e - 03 ± 2.46e - 04 ± 3.02e - 04 


1.70e-|-00 


2.00e-01 


1.84e - 03 ± 1.96e - 04 ± 1.47e - 04 



Table A. 3: 7r+ spectra in p-|-p coUisions. The unit of pt and pt width is GeV/c. 



Pt 


width 


M.B. 


0%-20% 


3.50e-01 


l.OOc-01 


3.20e -1- 00 ± 3.73e - 02 ± 2.56e - 01 


6.77e -1- 00 ± 9.88e - 02 ± 5.42e - 01 


4.50e-01 


l.OOc-01 


1.82e -1- 00 ± 2.46e - 02 ± 1.45e - 01 


3.92e -1- 00 ± 6.60e - 02 ± 3.13e - 01 


5.50e-01 


l.OOc-01 


1.08e -f 00 ± 1.72e - 02 ± 8.60e - 02 


2.42e -1- 00 ± 4.71e - 02 ± 1.94e - 01 


6.50e-01 


l.OOc-01 


6.70e - 01 ± 1.30e - 02 ± 5.36e - 02 


1.49e 00 ± 3.51e - 02 ± 1.20e - 01 


7.50e-01 


l.OOe-01 


4.05e - 01 it 9.06e - 03 ± 3.24e - 02 


9.22e - 01 ± 2.50e - 02 ± 7.38e - 02 


8.50e-01 


l.OOe-01 


2.59e - 01 ± 6.89e - 03 ± 2.07e - 02 


6.13e - 01 ± 1.96e - 02 ± 4.91e - 02 


9.50e-01 


l.OOe-01 


1.68e - 01 ± 5.38e - 03 ± 1.34e - 02 


3.85e - 01 ± 1.48e - 02 ± 3.08e - 02 


1.05e-|-00 


l.OOe-01 


1.16e - 01 ± 4.23e - 03 ± 9.29e - 03 


2.70e - 01 ± 1.18e - 02 ± 2.16e - 02 


1.15e-^00 


l.OOc-01 


7.63e - 02 ± 3.28e - 03 ± 6.11e - 03 


1.71e - 01 ± 8.83e - 03 ± 1.37e - 02 


1.30e+00 


2.00C-01 


4.36e - 02 ± 1.66e - 03 ± 3.49e - 03 


9.67e - 02 ± 4.45e - 03 ± 7.73e - 03 


1.50C+00 


2.00C-01 


2.04e - 02 ± 1.07e - 03 ± 1.63e - 03 


4.94e - 02 ± 3.06e - 03 ± 3.96e - 03 


1.70C+00 


2.00C-01 


1.03e - 02 ± 7.23e - 04 ± 8.26e - 04 


2.54e - 02 ± 2.07e - 03 ± 2.03e - 03 



Table A. 4: vr spectra in minimum-bias and 0-20% d-l-Au collisions. The unit of pr 
and Pt width is GeV/ c. 
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PT 


Pt width 


20%-40% 


40%-100% 


3.50e-01 


l.OOe-01 


4.81e + 00 ± 7.11e - 02 ± 3.85e - 01 


2.00e -1- 00 ± 2.95e - 02 ± 1.60e - 01 


4.50e-01 


l.OOe-01 


2.76e -1- 00 ± 4.73e - 02 ± 2.21e - 01 


1.12e -1- 00 ± 1.93e - 02 ± 8.95e - 02 


5.50C-01 


l.OOc-01 


1.65e + 00 ± 3.28e - 02 ± 1.32e - 01 


6.43e - 01 ± 1.31e - 02 ± 5.14e - 02 


6.50C-01 


l.OOc-01 


1.03e + 00 ± 2.47e - 02 ± 8.23e - 02 


4.04e - 01 ± 9.83e - 03 ± 3.23e - 02 


7.50C-01 


l.OOc-01 


6.15e - 01 ± 1.71e - 02 ± 4.92e - 02 


2.40e - 01 ± 6.80e - 03 ± 1.92e - 02 


8.50C-01 


l.OOe-01 


3.92e - 01 ± 1.30e - 02 ± 3.13e - 02 


1.51e - 01 it 5.10e - 03 ± 1.21e - 02 


9.50C-01 


l.OOc-01 


2.61e - 01 ± 1.03e - 02 ± 2.09e - 02 


9.62e - 02 it 3.91e - 03 it 7.70e - 03 


1.05e+00 


l.OOe-01 


1.81r - 01 ± 8.11k - 03 ± l.L^e - 02 


6.6.5(' - 02 ± 3.07(- - 03 ± .5.32c - 03 


l.l.jc+OO 


l.OOc-Ol 


1.20c - 01 ± ().2()c - 03 ^ <)..58(. - 03 


i.23c - 02 ^ 2.31( - 03 = 3.38c - 03 


1.30e+00 


2.00e-01 


6.82e - 02 ± 3.18e - 03 ± 5.46e - 03 


2.42e - 02 it 1.18e - 03 it 1.93e - 03 


1.50e+00 


2.00e-01 


3.24e - 02 ± 2.13e - 03 ± 2.59e - 03 


1.05e - 02 it 7.69e - 04 i 8.42e - 04 


1.70e+00 


2.00e-01 


1.55e - 02 ± 1.34e - 03 ± 1.24e - 03 


5.02e - 03 it 1.06e - 03 it 4.01e - 04 



Table A.5: 7r~ spectra in 20-40% and 40-100% d-^Au collisions. The unit of pr and 
Pt width is GeV/c. 



Pt 


Pt width 


P+P 


3.50e-01 


l.OOe-01 


9.71e - 01 it 1.22e - 02 ± 7.76e - 02 


4.50e-01 


l.OOe-01 


5.47e - 01 it 8.02e - 03 ± 4.37e - 02 


5.50e-01 


l.OOe-01 


3.09e - 01 it 5.38e - 03 it 2.47e - 02 


6.50e-01 


l.OOc-01 


1.84e - 01 ± 3.89e - 03 ± 1.47e - 02 


7.50C-01 


l.OOc-01 


l.OOe - 01 it 2.50e - 03 ± 8.01e - 03 


8.50e-01 


l.OOe-01 


6.36e - 02 it 1.89e - 03 ± 5.09e - 03 


9.50C-01 


l.OOc-01 


3.80e - 02 it 1.38e - 03 ± 3.04e - 03 


1.05e+00 


l.OOc-01 


2.44e - 02 it 1.03e - 03 ± 1.95e - 03 


1.15e-t-00 


l.OOe-01 


1.57e - 02 it 7.87e - 04 it 1.25e - 03 


1.30e-)-00 


2.00e-01 


8.70e - 03 ± 4.02e - 04 i 6.96e - 04 


1.50e-f00 


2.00e-01 


3.62e - 03 it 2.45e - 04 i 2.90e - 04 


1.70e-|-00 


2.00e-01 


1.69e - 03 it 1.76e - 04 it 1.35e - 04 



Table A. 6: tt spectra in p-|-p collisions. The unit of pt and pt width is GeV/c. 



Pt 


Pt width 


M.B. 


0%-20% 


4.57C-01 


l.OOe-01 


2.41e - 01 ± 9.47e - 03 ± 1.93e - 02 


4.83e - 01 it 2.75e - 02 it 3.86e - 02 


5.56C-01 


l.OOc-01 


1.87e - 01 ± 7.53e - 03 ± 1.49e - 02 


3.87e - 01 it 2.14e - 02 it 3.10e - 02 


6.55C-01 


l.OOe-01 


1.33e - 01 ± -S.ldc - 03 ± 1.07e - 02 


2.81e - 01 ± 1.53e - 02 ± 2.25e - 02 


7.54C-01 


l.OOe-01 


l.Ole - 01 ± 3.18e - 03 ± 8.07e - 03 


2. Me - 01 ± 1.09e - 02 it 1.72e - 02 


8.54C-01 


l.OOe-01 


6.97e - 02 ± 2.48e - 03 it 5.57e - 03 


1.58e - 01 ± 8.74e - 03 it 1.27e - 02 


9.54e-01 


l.OOe-01 


5.01e - 02 it 1.95e - 03 it 4.01e - 03 


9.95e - 02 it 6.30e - 03 it 7.96e - 03 


1.05e-|-00 


l.OOe-01 


3.77e - 02 it 1.65e - 03 it 3.02e - 03 


8.92e - 02 it 5.83e - 03 it 7.14e - 03 


1.15e-f00 


l.OOe-01 


2.73e - 02 it 1.31e - 03 it 2.18e - 03 


5.61e - 02 it 4.30e - 03 it 4.49e - 03 


1.30e-|-00 


2.00e-01 


1.78e - 02 it 7.25e - 04 i 1.42e - 03 


4.00e - 02 it 2.49e - 03 it 3.20e - 03 


1.50e-f00 


2.00C-01 


9.12e - 03 it 5.25e - 04 it 7.30e - 04 


1.95e - 02 it 1.89e - 03 it 1.56e - 03 


1.70e+00 


2.00e-01 


4.68e - 03 it 3.96e - 04 i 3.75e - 04 





Table A. 7: spectra in minimum-bias and 0-20% d-|-Au collisions. The unit of Pt 
and Pt width is GeV/ c. 
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Pt 


Pt width 


20%-40% 


40%- 100% 


4.57e-01 


l.OOe-01 


3.63e - 01 ± 2.06e - 02 ± 2.91e - 02 


1.27e - 01 ± 7.70e - 03 ± 1.02e - 02 


5.56e-01 


l.OOe-01 


2.74e - 01 ± 1.54e - 02 ± 2.19e - 02 


1.07e - 01 ± 6.16e - 03 ± 8.60e - 03 


6.55e-01 


l.OOe-01 


2.01e - 01 ± l.lle - 02 ± 1.60e - 02 


7.37e - 02 ± 4.26e - 03 ± 5.89e - 03 


7.54e-01 


l.OOe-01 


1.58e - 01 ± 8.18e - 03 ± 1.27e - 02 


5.25e - 02 ± 2.99e - 03 ± 4.20e - 03 


8.54e-01 


l.OOe-01 


l.Ole - 01 ± 6.16e - 03 ± 8.06e - 03 


3.73e - 02 ± 2.35e - 03 ± 2.98e - 03 


9.54C-01 


l.OOc-01 


7.97e - 02 ± 4.89e - 03 ± 6.37e - 03 


2.80e - 02 ± 1.86e - 03 ± 2.24e - 03 


1.05C+00 


l.OOc-01 


5.34e - 02 ± 3.83e - 03 ± 4.27e - 03 


1.91e - 02 ± 1.47e - 03 ± 1.52e - 03 


1.15C+00 


l.OOc-01 


3.91e - 02 ± 3.12e - 03 ± 3.13e - 03 


1.36e - 02 ± l.lSe - 03 ± 1.09e - 03 


1.30e+00 


2.00e-01 


2.51e - 02 ± 1.67e - 03 ± 2.01e - 03 


8.65e - 03 ± 6.83e - 04 ± 6.92e - 04 


1.50e+00 


2.00e-01 


1.32e - 02 ± 1.19e - 03 ± 1.05e - 03 


4.69e - 03 ± 5.07e - 04 ± 3.75e - 04 



Tabic A.8: K+ spectra in 20-40% and 40-100% d+Au collisions. The unit oipr and 
Pt width is GeV/c. 



Pt 


Pt width 


P+P 


4.57e-01 


l.OOe-01 


6.47e - 02 ± 3.01e - 03 ± 5.18e - 03 


5.56e-01 


l.OOe-01 


4.32e - 02 ± 2.10e - 03 ± 3.45e - 03 


6.55e-01 


l.OOe-01 


3.18e - 02 ± 1.51e - 03 ± 2.54e - 03 


7.54C-01 


l.OOc-01 


2.18e - 02 ± 9.70e - 04 ± 1.74e - 03 


8.54C-01 


l.OOc-01 


1.57e - 02 ± 7.69e - 04 ± 1.25e - 03 


9.54C-01 


l.OOc-01 


9.92e - 03 ± 5.60e - 04 ± 7.93e - 04 


1.05C+00 


l.OOc-01 


7.17e - 03 ± 4.74e - 04 ± 5.74e - 04 


1.15C+00 


l.OOo-Ol 


5.60e - 03 ± 4.18e - 04 ± 4.48e - 04 


1.30e+00 


2.00e-01 


3.73e - 03 ± 2.72e - 04 ± 2.98e - 04 


1.50e-)-00 


2.00e-01 


1.81e - 03 ± 1.18e - 04 ± 1.45e - 04 



Table A. 9: spectra in p-hp collisions. The unit of pr and pr width is GeV/c. 



PT 


Pt width 


M.B. 


0%-20% 


4.57C-01 


l.OOc-01 


2.34e - 01 ± 9.47e - 03 ± 1.87e - 02 


4.87e - 01 ± 2.81e - 02 ± 3.90e - 02 


5.56C-01 


l.OOc-01 


1.92e - 01 ± 7.98e - 03 ± 1.54e - 02 


4.01e - 01 ± 2.28e - 02 ± 3.21e - 02 


6.55e-01 


l.OOe-01 


1.39e - 01 ± 5.81e - 03 ± 1.12e - 02 


2.91e - 01 ± 1.69e - 02 ± 2.33e - 02 


7.54e-01 


l.OOe-01 


9.21e - 02 ± 3.10e - 03 ± 7.37e - 03 


1.90e - 01 ± 1.04e - 02 ± 1.52e - 02 


8.54e-01 


l.OOe-01 


7.11e - 02 ± 2.61e - 03 ± 5.69e - 03 


1.45e - 01 ± 8.76e - 03 ± 1.16e - 02 


9.54e-01 


l.OOe-01 


4.70e - 02 ± 1.94e - 03 ± 3.76e - 03 


9.49e - 02 ± 6.34e - 03 ± 7.59e - 03 


1.05e+00 


l.OOc-01 


3.11e - 02 ± 1.49e - 03 ± 2.48e - 03 


6.56e - 02 ± 5.07e - 03 ± 5.25e - 03 


1.15e+00 


l.OOc-01 


2.28e - 02 ± 1.21e - 03 ± 1.82e - 03 


4.74e - 02 ± 4.00e - 03 ± 3.79e - 03 


1.30C+00 


2.00C-01 


1.61e - 02 ± 7.10e - 04 ± 1.29e - 03 


3.70e - 02 ± 2.41e - 03 ± 2.96e - 03 


1.50e+00 


2.00e-01 


9.47e - 03 ± 5.54e - 04 ± 7.58e - 04 


2.01e - 02 ± 1.78e - 03 ± 1.61e - 03 


1.70e-|-00 


2.00e-01 


4.44e - 03 ± 4.19e - 04 ± 3.55e - 04 





Table A. 10: K spectra in minimum- bias and 0-20% d-l-Au collisions. The unit of 
Pt and pt width is GeV/ c. 
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PT 


Pt width 


20%-40% 


40%-100% 


4.57e-01 


l.OOe-01 


3.24e - 01 ± 1.95e - 02 ± 2.59e - 02 


1.39e - 01 ± 8.25e - 03 ± l.lle - 02 


5.56e-01 


l.OOe-01 


2.66e - 01 ± 1.58e - 02 ± 2.13e - 02 


1.13e - 01 ± 6.63e - 03 ± 9.03e - 03 


6.55C-01 


l.OOc-01 


1.99e - 01 ± l.lSe - 02 ± 1.59e - 02 


7.50e - 02 ± 4.61e - 03 ± 6.00e - 03 


7.54C-01 


l.OOc-01 


1.36e ~ 01 ± 7.60e - 03 ± 1.08e - 02 


5.62e - 02 ± 3.16e - 03 ± 4.50e - 03 


8.54C-01 


l.OOc-01 


1.13e - 01 ± 6.64e - 03 ± 9.04e - 03 


3.78e - 02 ± 2.47e - 03 ± 3.03e - 03 


9.54C-01 


l.OOe-01 


7.35e - 02 ± 4.83e - 03 ± 5.88e - 03 


2.55e - 02 ± 1.82e - 03 ± 2.04e - 03 


1.05C+00 


l.OOc-01 


4.51e - 02 ± 3.57e - 03 ± 3.61e - 03 


1.93e - 02 it 1.51e - 03 ± 1.54e - 03 


1.15e+00 


l.OOe-01 


3.03r - 02 ± 2.71k - 03 ± 2.13(- - 03 


1.36c - 02 ± 1.20(- - 03 ± 1.09c - 03 


1.30C+00 


2.()0c-01 


2.:5.'')C - 02 ± l.(i()C - 03 ^ 1.88(. - 03 


7.8 Ic - 03 ^ (i. lit - 01 = (>.27(; - 04 


1.50e+00 


2.00e-01 


1.05e - 02 ± 1.20e - 03 ± 8.36e - 04 


5.10e - 03 ± 5.72e - 04 ± 4.08e - 04 



Table A.ll: K' spectra in 20-40% and 40-100% d+Au collisions. The unit of pt and 
Pt width is GeV/c. 



PT 


Pt width 


P+P 


4.57e-01 


l.OOe-01 


5.99e - 02 ± 2.91e - 03 ± 4.79e - 03 


5.56e-01 


l.OOe-01 


4.79e - 02 ± 2.35e - 03 ± 3.83e - 03 


6.55e-01 


l.OOe-01 


3.06e - 02 ± 1.58e - 03 ± 2.45e - 03 


7.54e-01 


l.OOe-01 


2.15e - 02 ± 9.81e - 04 ± 1.72e - 03 


8.54e-01 


l.OOe-01 


1.46e - 02 ± 7.67e - 04 ± 1.17e - 03 


9.54e-01 


l.OOc-01 


l.Ole - 02 ± 5.86e - 04 ± 8.09e - 04 


1.05e+00 


l.OOc-01 


7.87e - 03 ± 5.22e - 04 ± 6.30e - 04 


1.15e+00 


l.OOc-01 


5.25e - 03 ± 4.27e - 04 ± 4.20e - 04 


1.30e-t-00 


2.00e-01 


3.42e - 03 ± 2.57e - 04 ± 2.74e - 04 


1.50e-|-00 


2.00e-01 


1.75e - 03 ± 1.61e - 04 ± 1.40e - 04 



Table A. 12: K spectra in p+p collisions. The unit of pr and pt width is GeV/c. 



PT 


PT width 


M.B. 


0%-20% 


4.68e-01 


l.OOe-01 


1.88e - 01 ± 1.87e - 02 ± 2.44e - 02 


3.96e - 01 ± 4.51e - 02 ± 5.15e - 02 


5.63e-01 


l.OOe-01 


1.35e - 01 ± l.lle - 02 ± 1.75e - 02 


2.80e - 01 ± 2.67e - 02 ± 3.64e - 02 


6.61e-01 


l.OOe-01 


1.07e - 01 ± 8.89e - 03 ± 1.39e - 02 


2.10e - 01 ± 2.01e - 02 ± 2.73e - 02 


7.59e-01 


l.OOe-01 


8.02e - 02 ± 6.78e - 03 ± 1.04e - 02 


1.74e - 01 ± 1.67e - 02 ± 2.27e - 02 


8.58C-01 


l.OOc-01 


5.82e - 02 ± 5.20e - 03 ± 7.57e - 03 


1.25e - 01 ± 1.27e - 02 ± 1.63e - 02 


9.57C-01 


l.OOc-01 


4.45e - 02 ± 4.85e - 03 ± 5.78e - 03 


1.05e - 01 ± 1.25e - 02 ± 1.37e - 02 


l.OOc+OO 


l.OOc-01 


3.63e - 02 ± 4. lie - 03 ± 2.90e - 03 


8.16e - 02 ± 1.02e - 02 ± 6.53e - 03 


1.16C+00 


l.OOe-01 


2.82e - 02 ± 2. lie - 03 ± 2.26e - 03 


5.08e - 02 ± 5.11e - 03 ± 4.06e - 03 


1.31e+00 


2.00e-01 


1.86e - 02 ± l.Me - 03 ± 1.49e - 03 


4.14e - 02 ± 3.21e - 03 ± 3.31e - 03 


1.50e-l-00 


2.00e-01 


1.02e - 02 ± 7.85e - 04 ± 8.14e - 04 


2.46e - 02 ± 2.34e - 03 ± 1.97e - 03 


1.70e-|-00 


2.00e-01 


5.64e - 03 ± 3.24e - 04 ± 4.51e - 04 


1.21e - 02 ± 1.03e - 03 ± 9.66e - 04 


1.90e-f00 


2.00e-01 


3.14e - 03 ± 2.33e - 04 ± 2.51e - 04 


8.08e - 03 ± 8.38e - 04 ± 6.46e - 04 


2.25e-|-00 


5.00e-01 


1.39e - 03 ± 9.47e - 05 ± 1.12e - 04 


3.37e - 03 ± 3.21e - 04 ± 2.70e - 04 


2.75e+00 


S.OOc-Ol 


2.75e - 04 ± 3.88e - 05 ± 2.20e - 05 


6.19e - 04 ± 1.26e - 04 ± 4.95e - 05 


3.50e+00 


l.OOc+00 


8.13e - 05 ± 1.37e - 05 ± 6.50e - 06 





Table A. 13: p spectra in minimum-bias and 0-20% d-|-Au collisions. The unit of Pt 
and Pt width is GeV/ c. 
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Vt 


Pt width 


20%-40% 


40%-100% 


4.68e-01 


l.OOe-01 


2.78e - 01 it 3.24e - 02 ± 3.61e - 02 


1.21e - 01 ± 1.40e - 02 ± 1.57e - 02 


5.63e-01 


l.OOe-01 


1.90e - 01 ± 1.84e - 02 ± 2.48e - 02 


7.84e - 02 ± 7.61e - 03 ± 1.02e - 02 


6.61e-01 


l.OOe-01 


1.56e - 01 ± 1.49e - 02 ± 2.03e - 02 


6.18e - 02 ± 5.96e - 03 ± 8.03e - 03 


7.59e-01 


l.OOe-01 


1.18e - 01 ± 1.14e - 02 ± 1.53e - 02 


4.84e - 02 ± 4.71e - 03 ± 6.29e - 03 


8.58e-01 


l.OOe-01 


8.54e - 02 ± 8.68e - 03 ± l.lle - 02 


2.98e - 02 ± 3.12e - 03 ± 3.87e - 03 


9.57e-01 


l.OOc-01 


6.45e - 02 ± 7.77e - 03 ± 8.38e - 03 


2.39e - 02 ± 2.93e - 03 ± 3.11e - 03 


1.06e+00 


l.OOc-01 


5.48e - 02 ± 6.88e - 03 ± 4.38e - 03 


2.12e - 02 ± 2.69e - 03 ± 1.69e - 03 


1.16C+00 


l.OOc-01 


4.21e - 02 ± 3.97e - 03 ± 3.37e - 03 


1.52e - 02 ± 1.49e - 03 ± 1.22e - 03 


l.Slc+OO 


2.00C-01 


2.68e - 02 ± 2.09e - 03 ± 2.15e - 03 


9.80e - 03 ± 7.91e - 04 ± 7.84e - 04 


1.50C+00 


2.00C-01 


1.41e - 02 ± 1.38e - 03 ± 1.13e - 03 


5.05e - 03 ± 5.17e - 04 ± 4.04e - 04 


i.roc+oo 


2.00C-01 


7.69e - 03 ± 7.02e - 04 ± 6.15e - 04 


2.54e - 03 ± 2.55e - 04 ± 2.03e - 04 


1.90e+00 


2.00e-01 


4.12e - 03 ± 5.31e - 04 ± 3.30e - 04 


1.43e - 03 ± 2.03e - 04 ± 1.15e - 04 


2.25e+00 


5.00e-01 


2.30(' - 03 ± 2.28k - 04 ± 1.84(- - 04 


6.31(' - 04 ± 7.57(- - 05 ± 5.05c - 05 


2.75('+00 


•l.OOc-Ol 


i.3()(. - 1 ± 1.13c -01 = 3. 19(. - 0.5 


l.()2t - 01 = 1.1 k - 05 = 1.29c - 05 



Table A. 14: p spectra in 20-40% and 40-100% d-l-Au collisions. The unit of pr and 
Pt width is GeV/ c. 



PT 


Pt width 


P+P 


4.68e-01 


l.OOe-01 


4.51e - 02 ± 5.24e - 03 ± 5.86e - 03 


5.63e-01 


l.OOe-01 


3.33e - 02 ± 3.29e - 03 ± 4.33e - 03 


6.61e-01 


l.OOe-01 


2.60e - 02 ± 2.57e - 03 ± 3.38e - 03 


7.59e-01 


l.OOe-01 


2.03e - 02 ± 2.12e - 03 ± 2.64e - 03 


8.58e-01 


l.OOe-01 


1.14e - 02 ± 1.30e - 03 ± 1.49e - 03 


9.57C-01 


l.OOc-01 


8.93e - 03 ± 1.20e - 03 ± 1.16e - 03 


1.06C+00 


l.OOc-01 


6.98e - 03 ± l.OOe - 03 ± 5.59e - 04 


1.16C+00 


l.OOc-01 


4.68e - 03 ± 5.95e - 04 ± 3.75e - 04 


1.31C+00 


2.00C-01 


2.90e - 03 ± 3.14e - 04 ± 2.32e - 04 


1.50C+00 


2.00C-01 


1.34e - 03 ± 2.04e - 04 ± 1.07e - 04 


1.70e+00 


2.00e-01 


6.61e - 04 ± 6.44e - 05 ± 5.29e - 05 


1.90e-|-00 


2.00e-01 


4.73e - 04 ± 5.58e - 05 ± 3.78e - 05 


2.25e-|-00 


5.00e-01 


1.84e - 04 ± 2.23e - 05 ± 1.47e - 05 


2.75e-|-00 


5.00e-01 


3.06e - 05 ± 7.26e - 06 ± 2.45e - 06 



Table A. 15: p spectra in p-|-p collisions. The unit of pt and px width is GeV/c. 



PT 


Pt width 


M.B. 


0%-20% 


4.68C-01 


l.OOc-01 


1.22e - 01 it 7.62e - 03 ± 1.59e - 02 


2.57e - 01 it 2.17e - 02 it 3.34e - 02 


5.63C-01 


l.OOc-01 


1.08e - 01 ± 5.87e - 03 ± 1.40e - 02 


2.26e - 01 ± 1.62e - 02 ± 2.94e - 02 


6.61C-01 


l.OOc-01 


8.86e - 02 ± 4.69e - 03 ± 1.15e - 02 


1.73e - 01 ± 1.23e - 02 ± 2.25e - 02 


7.59C-01 


l.OOc-01 


6.40e - 02 ± 3.47e - 03 it 8.32e - 03 


1.39e - 01 ± 9.76e - 03 it 1.81e - 02 


8.58e-01 


l.OOe-01 


5.51e - 02 it 3.04e - 03 it 7.16e - 03 


1.19e - 01 it 8.58e - 03 it 1.54e - 02 


9.57e-01 


l.OOe-01 


3.80e - 02 ± 2.30e - 03 ± 4.94e - 03 


8.99e - 02 ± 6.92e - 03 it 1.17e - 02 


1.06e-|-00 


l.OOe-01 


3.27e - 02 ± 1.41e - 03 it 2.62e - 03 


7.36e - 02 it 5.01e - 03 it 5.89e - 03 


1.16e-|-00 


l.OOe-01 


2.36e - 02 it 1.12e - 03 it 1.89e - 03 


4.24e - 02 it 3.50e - 03 it 3.39e - 03 


1.31e+00 


2.00C-01 


1.55e - 02 ± 6.16e - 04 i 1.24e - 03 


3.48e - 02 it 2.14e - 03 it 2.78e - 03 


1.50e+00 


2.00C-01 


8.00e - 03 it 4.04e - 04 i 6.40e - 04 


1.93e - 02 it 1.46e - 03 it 1.54e - 03 


1.70C+00 


2.00C-01 


4.78e - 03 ± 2.97e - 04 i 3.82e - 04 


l.Ole - 02 it 9.74e - 04 it 8. lie - 04 


1.90C+00 


2.00C-01 


2.56e - 03 ± 2.04e - 04 it 2.05e - 04 


4.99e - 03 it 6.54e - 04 it 3.99e - 04 


2.25e+00 


5.00e-01 


1.06e - 03 ± 7.97e - 05 ± 8.49e - 05 


2.69e - 03 it 2.85e - 04 it 2.15e - 04 


2.75e+00 


5.00e-01 


3.32e - 04 it 4.61e - 05 it 2.66e - 05 


5.86e - 04 it 1.35e - 04 i 4.69e - 05 


3.50e-|-00 


l.OOe+00 


7.89e - 05 it 1.53e - 05 it 6.31e - 06 





Table A. 16: p spectra in minimum-bias and 0-20% d-l-Au collisions. The unit of px 
and Pt width is GeV/ c. 
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PT 


Pt width 


20%-40% 


40%-100% 


4.68C-01 


l.OOc-01 


1.63e - 01 ± 1.45e - 02 ± 2.12e - 02 


7.50e - 02 ± 6.44e - 03 ± 9.75e - 03 


5.63C-01 


l.OOc-01 


1.65e - 01 ± 1.20e - 02 ± 2.14e - 02 


6.04e - 02 ± 4.55e - 03 ± 7.85e - 03 


6.61C-01 


l.OOc-01 


1.24e - 01 ± 8.94e - 03 ± 1.62e - 02 


5.51e - 02 ± 3.89e - 03 ± 7.16e - 03 


7.59e-01 


l.OOe-01 


9.35e - 02 ± 6.78e - 03 ± 1.22e - 02 


3.61e - 02 ± 2.68e - 03 ± 4.69e - 03 


8.58e-01 


l.OOe-01 


7.98e - 02 ± 5.96e - 03 ± 1.04e - 02 


2.85e - 02 ± 2.22e - 03 ± 3.70e - 03 


9.57e-01 


l.OOe-01 


5.22e - 02 ± 4.33e - 03 ± 6.78e - 03 


1.95e - 02 ± 1.68e - 03 ± 2.54e - 03 


1.06e+00 


l.OOe-01 


4.14e - 02 ± 3.20e - 03 ± 3.31e - 03 


1.82e - 02 ± 1.37e - 03 ± 1.45e - 03 


1.16e+00 


l.OOe-01 


3.61e - 02 ± 2.82e - 03 ± 2.89e - 03 


1.22e - 02 ± 1.04e - 03 ± 9.73e - 04 


1.31e+00 


2.00e-01 


2.32e - 02 ± 1.50e - 03 ± 1.85e - 03 


7.64e - 03 ± 5.47e - 04 ± 6.11e - 04 


1.50C+00 


2.00C-01 


1.13e - 02 ± 9.51e - 04 ± 9.04e - 04 


3.51e - 03 ± 3.35e - 04 ± 2.81e - 04 


1.70e+00 


2.00e-01 


6.41e - 03 ± 6.63e - 04 ± 5.13e - 04 


2.21e - 03 ± 2.48e - 04 ± 1.77e - 04 


1.90e+00 


2.00e-01 


4.64e - 03 ± 5.37e - 04 ± 3.71e - 04 


9.80e - 04 ± 1.70e - 04 ± 7.84e - 05 


2.25e+00 


5.00e-01 


1.46e - 03 ± 1.78e - 04 ± 1.17e - 04 


5.66e - 04 ± 7.10e - 05 ± 4.53e - 05 


2.75e+00 


5.00e-01 


4.43e - 04 ± 1.38e - 04 ± 3.54e - 05 


1.29e - 04 ± 4.62e - 05 ± 1.03e - 05 



Tabic A. 17: p spectra in 20-40% and 40-100% d-l-Au collisions. The unit of and 
Pt width is GeV/c. 



PT 


Pt width 


P+P 


4.68C-01 


l.OOc-01 


3.74e - 02 ± 2.60e - 03 ± 4.86e - 03 


5.63C-01 


l.OOc-01 


2.84e - 02 ± 1.77e - 03 ± 3.70e - 03 


6.61C-01 


l.OOc-01 


2.27e - 02 it 1.38e - 03 ± 2.96e - 03 


7.59C-01 


l.OOc-01 


1.32e - 02 ± 8.65e - 04 ± 1.71e - 03 


8.58e-01 


l.OOc-01 


1.05e - 02 ± 7.20e - 04 ± 1.37e - 03 


9.57e-01 


l.OOe-Ol 


7.17e - 03 ± 5.40e - 04 ± 9.32e - 04 


1.06e-f00 


l.OOe-01 


5.60e - 03 ± 3.77e - 04 ± 4.48e - 04 


1.16e-|-00 


l.OOe-Ol 


3.74e - 03 ± 2.86e - 04 ± 2.99e - 04 


1.31e-|-00 


2.00e-01 


2.31e - 03 ± 1.49e - 04 ± 1.84e - 04 


1.50e+00 


2.00C-01 


9.69e - 04 ± 8.84e - 05 ± 7.75e - 05 


1.70e+00 


2.00C-01 


5.95e - 04 ± 6.65e - 05 ± 4.76e - 05 


1.90e-^00 


2.00C-01 


3.57e - 04 ± 4.79e - 05 ± 2.86e - 05 


2.25e-|-00 


5.00e-01 


1.12e - 04 ± 1.90e - 05 ± 8.93e - 06 


2.75e-)-00 


5.00e-01 


3.87e - 05 ± 1.07e - 05 ± 3.10e - 06 



Table A. 18: p spectra in p-|-p coUisions. The unit of pr and pr width is GeV/c. 



Appendix B 

How to make MRPC 



This appendix is based on the procedure of the MRPC production in USTC. I will 
introduce the material preparations and then the chamber installation. 

B.l Preparations 
B.1.1 Glass 

(1) Check the glass very carefully by eye. The glass with scrapes is not accepted. 

(2) Measure the size of the glass with the digital vernier cahper. The errors of 
the length and width are required to be within 0.1 mm. Measure the thickness in 
several different places. The precision of the thickness is required to be 0.01 mm 
for each glass. (3) Use the micrometer to measure the flatness, which is required to 
be less than 0.01 mm. Use the mirror and observe the stripes of interference. (4) 
Grind the edge and the corner of the glass, and clean it. The size of outer glass is 
78{width) x206{length) x l.l{height) mm? and the size of inner glass is 61x200x0.54 
mm? . 

B.l. 2 Graphite Layer 

(1) Stick the layer in the middle of the outer glass. Squeeze the air out. (2) Stick a 
small copper tape, which is for high voltage (HV) applying, on to the graphite layer. 
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which is in the middle of the long side, and 0~0.5 mm away from the edge of the 
glass. The size of the graphite layer is 74 x 202 mm^. The size of the copper tape 
(the rectangle with the round angle) is 6 x 10 mrn^. 

B.1.3 Mylar layer 

Cut the mylar layer, and see if there is any tiny holes or scrapes. If yes, don't use it. 
The size of mylar is 84 x 212 x 0.35 mm^. 

B.1.4 Honeycomb board 

Measure the size and flatness. The error of the length and width is required to be 
within 0.2 mm, the error of the thickness is required to be within 0.05 mm. The 
flatness is required to be within 0.1 mm. The size of honeycomb board is 84 x 208 x 4 
mm^. 

B.1.5 The printed circuit board (PCB) 

(1) Check the surface of the metal which is used as read-out strips carefully, and see 
the position of the HV-holes is right or not. Check the size of the metal holes, whose 
diameters are required to be larger than 0.9 mm. (2) Use double side tape to stick 
the PCB board with the Mylar. The size of the double side tape is the same as the 
PCB board. The length of the mylar is 1 mm longer than that of the PCB board. 
(3) Use sealing ion to open a 3 mm hole, the center of the hole is in the middle of 
the HV holes. The size of PCB is 94 x 210 x 1.5 mm^. The size of metal holes are 
1 mm. 

B.1.6 Lucite cylinder 

Use digital vernier caliper to measure the length of the Lucite cylinder. Clean it and 
stick a double side tape on one side. The size is 3 mm and 3.87 < length < 3.93 
mm. 
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B.1.7 Other stuff 

Besides, we also need pins, fish line and little plastic cannula. The size of the pin is 
2-2.1 cm long. The fish line is (f) 0.22 mm. The plastic cannula is 1 mm. One kind 
of the cannula is 7 mm long, and the other is 5.6 mm long. Table IkTI lists the main 
materials for 1 MRPC 

B.2 Installation 

B.2.1 The outer glass and mylar and PCB 

(1) stick the outer glass on to the center of the mylar. (2) Use the sealing ion to 
connect the HV conductive line with the copper tape. Apply the HV to measure the 
noise rate and dark current. (3) Between the mylar and outer glass edge, on each 
side, use silica gel to seal. Attention: keep the surface clean and smooth. Attention: 
If one side is done, wait till the silica gel becomes solid. (4) Stick pins. Seal the pins 
which are used for the fish line coiling, into the metal holes of the PCB board. (5) 
Use the inner glass to fix on the position of Lucite cylinders, and keep them away 
from the pins for fish line. Then stick the 8~10 Lucite cylinders onto the outer glass. 

B.2. 2 Inner glass and fish-line coiling 

(1) Pre-install. Don't use fish-line. Pay attention to adjust the position of the pins. 

(2) This is now the real installation and fish line coiling. Clean the outer glass and 
inner glass carefully, coil a loop of fish line, add a piece of glass, then coil another loop 
of fish line, add another piece of glass, and so on and so forth. Attention: clean the 
fish line before it coils, and blow the surface of glass to protect it from the dirt with 
nitrogen jet. (3) Another time for pre-installation. Pay attention to the position of 
the upper and lower electrodes and adjust the position of pins. (4) Paste 3140 RTV 
coating onto the surface of Lucite cylinders. (5) Connect the two electrodes. Make 
sure all the pins connect right into the metal holes. Then lay the whole flat, and put 
on a block which is 4 kilogram weight. (6) After 2 hours, stick the honeycomb. (7) 
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iiiatciial 


type 


character 


number 


source 


outer glass 


window glass 


1.1 mm thick, VR: 8.7 X 10^-^ ohm. cm 


2 


Shanghai 


inner glass 


window glass 


0.54 mm thick, VR: 8.5 X 10^^ ohm.cm 


5 


USA 


graphite layer 


T9149 


0.13 mm thick, SR ; 2M ohm/square 


2 


Japan 


Mylar 


MO 


0.35 mm thick 


2 


Dupont Corp. 


honeycomb 




4 mm thick 


2 


Shanghai 


PCB 


gold 


1.58 mm thick 


2 


Shenzhen 


copper tape 




0.08 mm thick 


2 


3M Comp. 


LC 


Lucite 


(f) 3 mm, 3. 9 mm long 


8-10 


processing 


pins (single) 


metal pin 


21.5 mm long 


14 




pins (pair) 


metal pin 


21.5 mm long 


12 




cannula 


F-plastic 


(f> 1.4 mm 


38 




fish line 


top line 


<l> 0.22 mm 




Switzerland 


DST 


9690 


0.13 mm thick 




3M Comp. 


silica gel 


CAF4 


high-voltage insulation 




Switzerland 



Table B.l: The material for 1 MRPC model. VR is the volume resistivity and SR is 
the surface resistivity. LC is the Lucite cyhnder. DST is the double side tape. 

Measure the thickness of the whole. Make sure the precision is within 0.05 mm. (8) 
Connect the conductive-line for the read-out strips, and then put the whole into a 
bag. Attention: the conductive-line should not be broken. 
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